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THE SPECTRUM AND COMPOSITION OF THE 
GASEOUS NEBULAE 
By I. S. BOWEN 


ABSTRACT 


I. The primary source of nebular luminosity is the emission of radiation on the re- 
turn of electrons after ionization by extreme ultra-violet starlight. Additional luminos- 
ity is produced by the following secondary mechanisms: 

1. Lonization by secondary radiation.—Additional atoms are ionized by extreme ultra- 
violet quanta emitted by the primary mechanism. 

2. Excitation by secondary radiation.—Line absorption of secondary quanta may ex- 
cite an atom of the same kind as the primary emitting atom, or any other atom that has 
absorption lines which coincide with lines of the emitting atom. This mechanism, made 
possible by certain chance coincidences between lines of //e 11, O m1, and N 11, explains 
the presence with a high intensity of a selected group of O m1 and N 11 lines while all 
other lines of these elements are missing. 

3. Excitation by secondary electrons ejected by the primary mechanism.—This is prac- 
tically the sole source of excitation of the forbidden lines. 

II. A study of nebular line intensities in the light of the foregoing processes indi- 
cates that // is the most abundant element and //e is the second. V, O, Ne, S—and 
possibly C and A—are present but are very much rarer. The lines of these heavier ele- 
ments are strong, not because the elements are very abundant but because they are 
able to make use of large sources of energy that are not available to the predominent 
H and He. Lines of F, Na, Si, P, Cl, K, and Ca are missing. 


I. MECHANISM OF EXCITATION 


It is now generally recognized that the emission spectrum of the 
diffuse and planetary nebulae is produced by energy received from 
an adjacent star. The idea that the luminosity of the nebulae is ex- 
cited by adjacent stars is one of some antiquity, appearing first to 
have been suggested by Sir William Herschel’ in his discussion of the 
so-called “nebulous stars” discovered by himself. On the basis of 


' Abridged Phil. Trans., 17, 25, 1791. 
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information available to him, he was, however, disposed to favor a 
conception of the relationship of star to nebula that was unfavorable 
to this idea. As a consequence, this viewpoint received little atten- 
tion for a century or more. Its revival in recent years has resulted 
from the accumulation of observational data and from develop- 
ments in the theory of radiation. The favorable consideration it is 
now accorded as a working hypothesis is especially due to Hubble’s 
comprehensive discussion of observational data provided largely by 
V. M. Slipher, W. H. Wright, H. D. Curtis, and himself.’ 

The primary mechanism causing this line emission is the ioniza- 
tion of the atoms of the nebula by the absorption of stellar radiation 
of a wave-length shorter, i.e., of a frequency greater, than that cor- 
responding to the ionization potential. Thus, an H atom can be 
ionized by all wave-lengths less than 912A, He 1 by those less than 
504 A, and He 11 by those less than 228 A, etc. The observed spectra 
of these atoms are then emitted on the return of the ejected electron. 
This mechanism was first suggested by D. H. Menzel’ and by H. 
Zanstra,’ the latter of whom has discussed quantitatively the rela- 
tive intensities of the star and of the nebula excited by it on the 
basis of this mechanism. 

In addition to this primary process, however, there is in operation 
a series of secondary mechanisms which materially change not only 
the total intensity of the nebular radiation but the type of spectrum 
emitted as well. In order to simplify the discussion of these secondary 
processes they will be considered in terms of the train of events 
which may follow the ionization by the primary mechanism of one 
type of atom, namely, He. At some time after ionization the 
He** ion picks up an electron, usually in some outer orbit, from which 
the electron returns to its normal state by a series of transitions. 
The first transitions are normally between outer orbits and give rise 
to the observed He 11 nebular lines. In every case, however, the last 
transition made ends in the normal state and corresponds to the 
emission of one of the series of extreme ultra-violet lines at 304, 256, 
243 A, etc., depending on whether the final jump was made from the 

2 Ap. J., 56, 162, 400, 1922. 

3 Pub. A.S.P., 38, 295, 1926. 

+Ad. J .°GS; 50,1927; 25... Ap., a, 1, 1931. 
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2p, 3p, or 4p orbit. One or more of the following secondary pro- 
cesses may then take place: 

1. Ionization by secondary radiation.*—Since these extreme ultra- 
violet lines of He 11 are below the critical wave-length necessary to 
ionize both H and He1, they may ionize one of these atoms and 
thereby cause the emission of its observed spectrum. Indeed, if He 1 
is ionized first, the radiations given off by it may in turn ionize an H 
atom. 

2. Excitation by secondary radiation.—As the extreme ultra-violet 
lines of He 11 listed above are transitions to the normal state of the 
atom, they are very strongly absorbed by other He 11 atoms. When 
this type of absorption occurs, the He u atom is, of course, excited 
to the same upper state that the original emitting atom was in before 
the ultra-violet quantum was given off. In general, this He 1 atom 
thus excited can now re-emit the quantum absorbed, or if the ex- 
cited state is above the 2p state, it may emit first a line of long 
wave-length, then another of the extreme ultra-violet lines. As 
pointed out by Zanstra,° a series of absorptions and re-emissions of 
this type soon results in the transformation of all of these extreme 
ultra-violet quanta into quanta of 304 A, plus some of the He 11 lines 
of longer wave-length which can immediately escape from the nebula, 
since they are not absorbed by the normal He 11 atom. Since many 
of these lines of longer wave-length fall in the observable range, this 
mechanism adds to the intensity of the observed He 1 lines. 

As was recently pointed out,’ a few of the O m1 and JN 11 lines can 
be excited by these extreme ultra-violet He 11 lines by this same 
mechanism, which depends, in this case, on certain chance coinci- 
dences between the wave-lengths of He 11, O u1, and N 111 lines, as 
illustrated in the simplified Grotrian diagrams of Figure 1. From 
these diagrams it is at once evident that because of the near coinci- 
dence of the s’p*3P,—s’p3d3P, line of O11 at 303.799 A with the 
303.780 A resonance line of He 1, the O m1 atom can be excited to 
the s’*p3dP, state by the absorption of this very intense He 11 line. 
From this level the electron can either return directly to the s*p?3P 
states, thus re-emitting the line at 304 A, or it can jump first to a 

51. S. Bowen, Ap. J., 67, 14, 1928. 

® Thid., 65, 56, 1927. 71. S. Bowen, Pub. A.S.P., 46, 146, 1934. 
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term of the s*p3p configuration, next to a term of the s*p3s configura- 
tion, and then to the normal state. Obviously, since the quantum is 
not lost when the 304 A line is re-emitted, the eventual result of sev- 
eral absorptions and re-emissions of this line is nearly always the 
production of two of the observed long wave-length O u1 lines fol- 
lowed by a quantum of 374 A. 
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Nui On REP Heil 





3kD 74 1% 
4647 
464229745 | | | 
VMS / 
22D) b 
Ip ia \ iy 


rf / | | 
& || | | 





\ 374434) /374.942 374436\ 303777] /303673 303780) 
\ | \ 
\ / | 
\ | 
\ // \ | 
452 \ / / \ | 
oat // \ f 
\ i} \ 
\ / \ 
. \ 79 \ rac | 
2pP Mw% 2pP Nz 1 2 | se a 


Fic. 1.—Terms of Nim, O1n, and Het taking part in secondary mechanism No. 2 


A similar near coincidence between one of these 374 A lines of 
O11 and the s’*p’*P,,—s’3d’D., line of NV ur enables the observed 
N 111 lines to be excited by this same process. 

In general, the absorption coefficient of an atom for radiation co- 
inciding with one of its own natural frequencies emitted on the return 
of an electron to its normal state (secondary mechanism No. 2) is 
very much greater than the absorption coefficient for radiation hav- 
ing a wave-length below its ionization limit (primary mechanism or 
secondary mechanism No. 1). This fact enables the secondary mech- 
anism (No. 2) to play a very important réle in determining the ex- 
citation and therefore the appearance of the nebular spectra. 

In the first place, the dimensions of the nebula are of the order of 
magnitude of the mean free path of a quantum before the primary 
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absorption with ionization occurs. Consequently this much higher 
absorption coefficient of the abundant H 11 atoms for their own 304, 
250 A, etc., radiations prevents these radiations from escaping 
directly from the nebulae. Instead, these quanta must first be ab- 
sorbed and re-emitted many times, and if they eventually escape at 
all, it is only after a slow outward diffusion to the boundary of the 
nebula. This, therefore, gives the other secondary processes a much 
larger chance to make use of the 304 A quanta before they escape 
from the nebula. 

In the second place, it enables the majority of the 304 A quanta to 
excite the O11 lines and they in turn the N 11 lines before these 
quanta are absorbed by the abundant H and He 1 atoms by means of 
the ionization process. Because of this, the intensities of the O10 
and N 111 lines are proportional to the rate of formation of the 304 A 
quanta, but are, through wide limits, independent of the abundance 
of Out and N 1 even though these atoms are much rarer than 
those of H and He. 

That the observed O 111 and N 111 permitted nebular lines are al- 
most entirely excited by this mechanism is shown in Table I, which 
gives all of the O ur and N 11 permitted nebular lines that W. H. 
Wright* has observed in his very extensive studies. For comparison 
the stronger O 111? and N 11"? lines that appear in laboratory sources 
are also listed. In O 11 it is seen that all possible transitions, in the 
observable range, from the 3d’P, state, and all transitions from the 
three states (s*p3p%S, *P.3D,) to which an electron in the s’*p3d3P, is 
most likely to jump, are definitely present. On the other hand, all 
other O 111 lines, with one very faint exception, are completely miss- 
ing, in spite of the fact that in the laboratory all of the other lines 
listed are much stronger than some that are observed in the nebulae. 
The one very faint exception at 3118 A (s’*p3p3S —s’p3d3P,) is doubt- 
less also excited by the same mechanism, since a relative velocity of 
go km/sec., which may rarely occur, is enough to cause the excita- 
tion of the s’*p3dP, state by changing the wave-length of the He line 

8 Pub. Lick Obs., 13, 193, 1918; Lick Obs. Bull., 17, 1, 1934; Pub. A.S.P., 46, 142 and 
280, 1934. 

9A. Fowler, Proc. R. Soc., A, 117, 317, 1928. 


10 L, J. Freeman, ibid., 121, 318, 1928. 
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(304.780 A) to the wave-length (304.693 A) of the s*p*3P, —s’p3d5P, 
line of O 111. 

Similarly, with V 1 the lines that may be excited by this sec- 
ondary mechanism are present in the nebulae while all others are 
missing. Owing to the fivefold greater absorption coefficient of the 
374.434 A (s’p’P,,—s?3d?D.,) line than that of the 374.442 A 
(s*p’?P,, —s’3d?D,,) line, an unusually large number of atoms are ex- 


TABLE I 


PERMITTED LINES OF O III AND N III 


Int. \ Neb- Int. \ Neb- 
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cited to the s*3d?D., state, with the result that the 4641 A line and 
the 4642 A line, which are not resolved in the nebular spectrum, 
should be about seven times as strong as 4634 A. Wright" has noted 
just such a discrepancy between the laboratory and the nebular in- 
tensities. In case of high relative velocities between the O 11 and 
N ut, this discrepancy in intensity between 4641 and 4634 A may 
be reduced by the absorption of the 2p’P,;—3d?D,, line of NV ut at 
374.204 A, which differs by only 0.039 A, equivalent to a velocity of 
32 km/sec., from another component of the 374 A O ur multiplet. 

Further support for this mechanism of excitation of the O m1 and 


™ Pub. A.S.P., 67, 149, 1934. 
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N u1 lines is given by the fact that the sizes of the monochromatic 
images of the nebulae given by these lines are very nearly equal to 
those obtained with the He 11 lines but differ widely from the image 
sizes of the O 111 forbidden lines. Likewise, the variation in the in- 
tensities of these lines in going from nebulae to nebulae parallels that 
of the He 11 lines much more closely than it does that of the forbidden 
lines of O m1. 

This mechanism may also play quite a role in determining the ap- 
pearance of the spectra of novae, although due to the much wider 
lines present there a larger number of lines of O 11 and N m1 may be 
excited by it, thus making the evidence for the mechanism less 
striking. If this is the case, the concentration of the 452 A N 1m line, 
emitted at the end of the series of processes, may become very great 
in these somewhat denser objects. It may be this high radiation 
density of the resonance lines of V 11 that causes the nitrogen to 
appear to “blow up” at certain stages of the development of the 
novae, as indicated by the very sudden and marked widening of the 
nitrogen lines.” 

It is possible that the abnomally high intensity of the 4686 A 
(35,P,D—4S,P,D,F) line of He 11 in the nebulae may partially, at 
least, be due to the fact that the high intensity of 304 A keeps an 
appreciable number of atoms in the 2s’S state, which is probably 
metastable. From this state it may be raised to the 4p’P state by 
absorption of the still more intense 1216 A resonance line of H. A 
relative velocity of 120 km/sec. between the H and the He 11 atoms 
is required, but in view of the lightness of both atoms this may not 
be too rare an occurrence. 

3. Excitation by secondary electrons.'’—Since in general the quan- 
tum of energy absorbed by the atom in the primary process has 
more than enough energy to ionize the atom, the electron is ejected 
with a considerable velocity. Thus, as will be discussed in the next 
section, the only atoms abundant enough to play any large réle in 
the primary process are H, He1, and He 11, whose ionization poten- 
tials are 13.5, 24.5, and 54.1 electron volts, respectively. Since the 
primary radiation from the star is continuously distributed in wave- 

2W.H. Wright, Pub. Lick Obs., 14, 84, 1921. 

131. S. Bowen, Ap. J., 67, 13, 1928. 
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length, it is evident that the electrons will come off with energies up 
to at least the differences between these, i.e., 10-30 electron volts. 
The same is true of the electrons ejected from H and Het by the 
secondary radiations 304, 374, 452, and 585 A. The electrons with 
these energies may then excite other atoms by impact. 

This mechanism is particularly effective in exciting the forbidden 
lines, since these represent transitions from metastable states whose 
excitation potentials are very low. Thus the nebular pair at 4959 and 
5007 A has an excitation potential of only 2.50 electron volts, while 
those of other forbidden lines range from 1.83 to 5.33 electron volts. 
In comparison with these, the lowest excitation potentials of the 
very abundant H, He1, and He atoms are 10.2, 19.7, and 40.6 
electron volts, respectively. Consequently the very large source of 
energy in the form of electrons having energies of 10 or less electron 
volts is only available for the excitation of the forbidden lines. The 
intensity of these lines is therefore independent, over a wide range, 
of the abundance of the elements involved in comparison with the 
predominant H and He. 

This mechanism is practically the sole source of excitation of the 
forbidden lines since, as will be shown in the next section, the per- 
mitted lines of the elements concerned are almost entirely missing. 
These permitted lines should be emitted first if the electron arrived 
in the metastable state by returning after ionization by the primary 
mechanism. 

II. CHEMICAL COMPOSITION 

From the foregoing discussion of the mechanisms of excitation of 
the nebular lines it is evident that the strength of the forbidden lines 
or of the permitted O u1 and N 111 lines excited by line absorption is 
not necessarily an indication of the high relative abundance of the ele- 
ments emitting them. In order to make even a rough estimate of the 
relative abundance of various elements it is therefore necessary to 
make comparisons among permitted lines excited by the primary 
mechanism or among forbidden lines. 

Table II gives all of the permitted nebular lines, found in 
Wright’s" lists, that are excited chiefly by the primary mechanism. 
The intensities listed are Wright’s values obtained with the quartz 


™4 See n. 8, also Pub. A.S.P., 32, 63, 1920. 
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spectrograph in the region below 3850 A, and with the glass spectro- 
graph above this wave-length. If we neglect the 4686 A line of He uy, 
which may be abnormally strong owing to the secondary mechanism 
discussed above, we see that hydrogen is the predominant element 
in the nebulae. While considerably weaker than those of H, the 
Hei and He 11 lines are also definitely and consistently present. Be- 
sides those of H and He, however, only four lines that might be 
identified as permitted lines are found. Two of them, 3935 and 
4416 A, were marked by Wright as only suspected, since they were 
very faint and were observed in only one nebula each. In the case of 
3935 A, provisionally identified as Ca 11, the nebula is an abnormal 
one (BD 30°3639) and evidence discussed later indicates that Ca is 
not present in the more normal nebulae. The other two lines, 4267 
and 4650 A, were also very weak and were found in but four and 
two nebulae, respectively. Since in each of the four cases but one or 
two lines are attributed to a given ion, the identifications are all very 
doubtful. 

It is therefore evident that all elements except H and He are very 
rare in the nebulae. Consequently the only reason that certain lines 
of the heavier elements, such as the forbidden lines or the fluorescent 
lines of O 111 and N 101, are strong is that they can be excited by large 
sources of energy that are not available to the predominant H and 
He. 

It is also evident that the permitted lines of the heavier elements 
are either missing or so weak and uncertain that they cannot satis- 
factorily be used to estimate the relative abundance of these atoms. 
This forces one to make use of the intensities of the forbidden lines 
which, while they are not comparable with permitted lines, may be 
compared among themselves to estimate roughly the relative abund- 
ance of those heavier elements that are capable of emitting forbidden 
lines in the observable range. That the forbidden lines are much 
more sensitive for the detection of the rare elements is at once obvi- 
ous. Thus none of the O 111 permitted lines that are excited by the 
primary mechanism is strong enough to appear, while its forbidden 
lines at 4959 and 5007 A are by far the strongest lines in the whole 
spectrum of many nebulae. 

Table III lists the predicted positions of all forbidden lines that 
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PERMITTED LINES EXCITED BY PRIMARY MECHANISM 
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any of the elements up to Ca are likely to emit in the observable 
range between 2900 and 10,000 A. The uncertainty indicated in the 
fourth column is the maximum tolerance that is allowable in the pre- 
dicted positions. In general, uncertainties below 1o A are due to the 
large errors in wave-number units that are inherent in the extreme 
ultra-violet measurements on which these predictions depend. 
Larger uncertainties than this occur only in cases where direct 
measurements are not available and the predicted wave-lengths de- 
pend on interpolations or extrapolations from nearby atoms of the 
same structure. On account of a much lower excitation potential, 
transitions from the 'D or ?D states will, in any given atom, be much 
stronger than transitions from the 'S or ’P states. 

The nebular wave-lengths and intensities are again taken from 
Wright’s list, except that three values in the extreme red are ob- 
tained from P. W. Merrill.’ The references in the first and third 
columns indicate the spectroscopic data used in making the predic- 


5 Pub. A.S.P., 40, 254, 1928. 
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* B. Edlén, Nova Acta Reg. Soc. Sci. Upsaliensis, 9, No. 6, 1934. 
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tJ. C. McLennan and J. H. McLeod, Proc. R. Soc., A, 115, 515 (direct observation). 
§ J. J. Hopfield, Phys. Rev., 37, 160, 1931 (direct observation). 


| F. Becker and W. Grotrian, Ergebnisse der Exakten Naturwissenschaften, 7, 65, 1928. 


"1. S. Bowen, Ap. J., 67, 1, 1928. 
** Auroral green line. Wave-length measured by H. D. Babcock, ibid., 57, 200, 1923. Identification by 
J. C. McLennan, J. H. McLeod, and R. Ruedy, Phil. Mag., 6, 558, 1928. 
tt I. S. Bowen, Phys. Rev., 36, 600, 1930; F. Paschen, Naturwissenschaften, 18, 752, 1930. 
tt Measurements by Merrill, Pub. A.S.P., 40, 254, 1928. 


§§ I. S. Bowen, Phys. Rev., 45, 82, 1934. 
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© B. Edlén and P. Swings, C.R., 198, 1842, 1934. 

*** J.C. Boyce, Phys. Rev., 46, 378, 1934. 
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ttt P. Swings and B. Edlén, C.R., 198, 1748, 1934; I. S. Bowen, Pub. A.S.P., 46, 145, 1934. 

§$$§ J. Soderquist, Nova Acta Reg. Soc. Sci. Upsaliensis, 9, No. 7, 1934. 

|| A. Fowler, Proc. R. Soc., A, 123, 422, 1920. 

© C. C. Kiess, J. of Research, Bur. of Stand., 8, 393, 1032. 

*#4* J. E. Ruedy, Phys. Rev., 44, 757, 1933. Ruedy’s identification of s?p1*S —s?p3(?P)4s'P has been 
changed to the line at 1782.25 A to make the relationships of terms of the s?p* configuration agree better 
with other atoms of similar structure. 

tttt M. Gilles, Ann. d. Phys., 15, 267, 1931. 

ttitI. S. Bowen, Phys. Rev., 46, 377, 1934. 

§§$$$ I. S. Bowen, Nature, 123, 450, 1929. 
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tions. The notations on the observed wave-lengths refer to the arti- 
cle in which the identifications were originally made. 

From the table it is evident that V, O, Ne, and S are definitely 
present and that O is considerably more abundant than the others. 
Argon may be present, although there are serious inconsistencies in 
the nebular behavior of the A Iv pair. Consequently, further obser- 
vations in the extreme red are needed to determine whether the 
other component of the A 111 and A v pairs are definitely present or 
absent. The evidence is not conclusive in the case of CI, since the 
only line, 4621 A, in the observed range has, according to E. A. Co- 
don’s” calculations, less than 1/200 of the intensity of the 'D—'S 
line. It should therefore not be strong enough to be observed with 
any reasonable concentration of C. Likewise, the only Mg and Al 
lines in this range are weak transitions from the ’P levels of stages of 
ionization so high that their production would be doubtful even 
though these elements were present. On the other hand, the lines of 
F, Na, Si, P, Cl, K, and Ca appear to be rather definitely missing 
and the abundance of these elements must therefore be small. 

Unfortunately, the spectra of the astronomically common metals 
between Ca and Cu have not been analyzed to an extent that a sys- 
tematic study of their forbidden lines is possible. However, for the 
first two stages of ionization of these elements the excitation poten- 
tials of the permitted lines are not much greater than those of the 
forbidden lines. Consequently, if these elements are present in these 
lower stages of ionization we should expect the permitted lines to be 
observed as well as the forbidden lines. The failure to observe these 
lines therefore seems to eliminate the possibility of any very large 
concentration of these heavier atoms, at least in low stages of ioniza- 
tion. Indeed, less than a dozen nebular lines remain unidentified, 
and all of these are of such low intensity that they are barely ob- 
servable in from one to four of the brightest nebulae. Consequently 
one can, in general, assign an abundance lower than that of NV, Ne, or 
S to all unidentified elements which have, in the observed range, 
either permitted or forbidden lines of low excitation potential, even 
though the position of these lines cannot be fixed at the present 
time. 


16 Ap. J., 79, 217, 1934. 
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From this qualitative survey one may conclude, therefore, that 
like most astronomical bodies the nebulae are made up largely of 
hydrogen. Helium is the next most abundant element. In the third 
place, but very much rarer than these, one finds O, or possibly C if 
its 4267 and 4650 A lines are correctly identified. Still rarer, but 
definitely present, are NV, Ne, S, and possibly A. In comparison 
with the sun’? there seems to be a tendency to favor the non-metals. 
Thus if, as is true in the sun, Na, K, and Ca were about equally 
common with N and ten times as abundant as S, we should doubt- 
less find their lines in the nebular spectrum. 

Unfortunately, our present knowledge of such factors as ultra- 
violet absorption coefficients, internal nebular velocities, nebular 
densities, velocity distribution of nebular electrons, etc., is so meager 
that a quantitative discussion of the abundance of the elements as a 
function of the intensities of nebular lines is not feasible. 


In conclusion the writer wishes to express his indebtedness to 
Professor W. H. Wright and to Dr. J. C. Boyce for the use of un- 
published data, and for various suggestions concerning observational 
material on which this article is based. 

CALIFORNIA INSTITUTE OF TECHNOLOGY 


October 1934 
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LUNAR RADIATION AS RELATED TO PHASE 
By EDISON PETTIT* 


ABSTRACT 

Methods of observing —An image of the moon, 0.87 mm maximum diameter, was 
formed by a concave mirror of 6.4 cm aperture and 8.8 cm focal length alternately 
upon the two receivers of a vacuum thermocouple. The planetary heat was separated 
from the light by a microscope cover-glass. The apparatus was calibrated by a standard 
lamp at 71.3 m distance. An air-mass meter was devised to read sec Z directly. 

Reduction of the measurements.—The measurements were reduced to the zenith with 
the atmospheric extinction observed for the night or, when only one or two observations 
were made, by using the mean extinctions, 0.11 mag. for reflected light and 0.15 mag. 
for planetary heat. Small corrections were applied for the planetary heat transmitted, 
and light reflected, by the microscope cover-glass screen, and the measurements were 
then reduced to mean distance of sun and moon. The mean curves were corrected to no 
atmosphere for comparison with the theoretical formula. 

Radiometric magnitude-——The radiometric magnitude of planetary heat per square 
second of arc for temperatures between 140° and 500° K is given by the empirical equa- 
tion 

1545 


mi +Am/) = T —2.6. 


The radiometric magnitude of the planetary heat from the full moon under normal con- 
ditions is —14.8, and that of the reflected light, — 13.4. 

Comparison with other data—The measurements of planetary heat made by the 
Earl of Rosse are 10-20 per cent lower than the present values for most phases of the 
moon, and lower still for reflected light, although the uncertainty of the corrections 
which should be applied makes the comparison difficult. No great amount of asym- 
metry such as has been reported by others was observed in the phase-intensity curve. 
The theoretical curve is everywhere higher than the observed curve. 

Transmission of ozone for planetary heat-—The transmission of an equivalent of 2 cm 
of ozone measured for the differential radiation between boiling water and liquid air was 
0.92, which indicates that atmospheric ozone absorbs 4 per cent of the planetary heat 
from the moon. 


It is well known that the radiation from the moon consists of two 
parts—reflected sunlight and radiant heat emitted by its warmed 
surface—the latter called planetary heat. A study of the radiation 
from the subsolar point and other particular areas on the moon was 
given in a previous paper.’ It is now proposed to study the radiation 
from the moon as a whole. 

More than sixty-five years ago this problem was studied by the 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 504. ' 
' Pettit and Nicholson, Mt. W. Contr., No. 392; Ap. J., 71, 102, 1930. 
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Earl of Rosse,’ using a 3-foot speculum reflector and two large 
thermopiles at atmospheric pressure upon which the lunar image was 
alternately focused by two secondary concaves. The light was sepa- 
rated from the planetary heat by means of a plate of glass the thick- 
ness of which is not stated. Our increased knowledge of the trans- 
mission coefficients of the atmosphere and the instruments, and of 
the proper methods of reducing the measurements, made it desirable 
to repeat this work with a more sensitive arrangement and under 
better atmospheric conditions. 

The reflected light from the moon has been studied visually by 
several observers, and by J. Stebbins and F. C. Brown’ with a 
selenium cell. The latter observers used an instrument which was sen- 
sitive throughout the visible spectrum but which, for the most of their 
observations, showed some preference for the region around J 0.7 yu. 
Their measures were made by exposing the cell, contained in a 
blackened box, directly to moonlight and determining the distance 
at which an amyl acetate lamp gave the same deflection. More 
recently, G. Rougier has made measurements of moonlight with a 
photo-electric cell and amplifier.‘ 

No such simple method can be used to measure the lunar radiation 
with a thermocouple, but the same result may be obtained by using 
an optical system which forms an image of the entire moon upon a 
vacuum thermocouple receiver without the secondary mirrors used 
by Lord Rosse. With the aid of a microscope cover-glass' the re- 
flected light may be separated from the planetary heat with a mini- 
mum of loss, thus giving two sets of data, viz., the march with phase 
angle of the whole reflected solar radiation and of the planetary heat. 
Since we thus do with a smaller instrument what was done on the 
planet Mercury with the 1oo-inch telescope,’ the results should be 
comparable. 

2 Proc. R. Soc. Lon., 27, 436, 1869; 29, 9, 1870; Phil. Trans. R. Soc. Lon., 163, 587, 
1873. 

3 Ap. J., 26, 326, 1907. 


4 Ann. de l’Observatoire de Strasbourg, 2, 205, 1933. This paper gives a considerable 
history of the measurements of moonlight. See also L’Astronomie, 48, 220 and 281, 


1934. 
5 Pettit and Nicholson, Pubs. A.S.P., 35, 194, 1923; Pop. Astr., 31, 657, 1923, and 


33, 299, 1925. 
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THE APPARATUS 

Figure 1 shows a diagram of the apparatus used in the investiga- 
tion. A thermocouple, 7, consisting of two junctions with receivers 
I mm square, separated 2.5 mm, is placed at the focus of a 6.4-cm 
parabolic mirror, A, of 8.8 cm focal length. Since the maximum 
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Fic. 1.—Diagram of thermoelectric lunar radiometer. The air-mass meter is shown 


at the left. 


diameter of the image of the moon formed by this mirror is 0.87 
mm, the receiver measures all the lunar radiation thus reflected. 
The thermocouple, 7, is mounted in a small evacuated cell, 2 cm in 
diameter, connected with a large-capacity sphere, Q, and a calcium- 
filled tube used to assist in maintaining the vacuum. A rock-salt 
window, R, admits the radiation; and a white celluloid screen, S, 
cemented just behind the thermocouple, 7, in the cell, assists in 
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finding the lunar image. The plug, P, containing the lead wires of 
the thermocouple, and the tube, O, are in the same plane, and there- 
fore cast a common shadow on the mirror, A. A microscope cover- 
glass, § mm thick, is mounted on the arm, U, actuated by a plunger 
magnet, .V, in such a way that it may be placed in front of the ther- 
mocouple or removed by operating the switch, Y. The arm, U, is 
also in the plane of O and P and therefore in the common shadow. 
The mirror, A, has a hole, 13 mm diameter, in its center. By means 
of a single-lens eyepiece, /, the lunar image, 0.87 mm diameter, may 
be set on the thermocouple receivers. The instrument is mounted 
on a square aluminum plate, YZ, attached to a second plate by a 
fulcrum, F, about which it may rotate in its own plane. A small 
angular displacement of this plate about F is given by the motor- 
operated crank, Z, working in a slot in the plate. The motor, V, 
also operates the commutator, C, which permits the crank to make 
only half a revolution, after which the motor circuit is broken 
through a relay. Reversing the switch, W, causes the crank to make 
another half-revolution. Hence, by operating the switch, W, the 
image of the moon may be moved from one junction of the thermo- 
couple to the other. A guiding telescope, G, with one thread in right 
ascension and two in declination, serves to set the lunar image on 
the thermocouple receivers in daytime sky when working at crescent 
phases. 
THE AIR-MASS METER 

A device, K, which may be called an ‘‘air-mass meter,” is also a 
part of the apparatus. It consists of a plummet, O, with a right- 
angle arm, J, balanced by the mass, //, and fixed by the fulcrum, F, 
to the spindle, AB, which may rotate about its length. The gradu- 
ated edge of the spindle is the axis. The lower edge of the pointer, 
PJ, at all times lies in the plane of the horizon and in that of the 
vertical circle passing through the object. #P is 1 cm, and the scale 
edge is graduated in centimeters and millimeters. Since AB is set 
parallel to the optical axis of the reflector, A, and guiding telescope, 
G, the intersection of PJ with the scale AB gives the value of sec Z, 
or the air-mass. Although the readings cannot be made closer than 
about 0.05 air-mass, they are useful in checking the reduction for 
air-mass on a large number of observations. 
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- 











THERMOELECTRIC LUNAR RADIOMETER 


omg 
Or 
ary 





LUNAR RADIATION AS RELATED TO PHASE 21 


The apparatus is mounted on the declination axis of the 6-inch 
telescope, replacing most of the counterweights (Pl. I). The tele- 
scope tube, J, is cut longitudinally and splayed in the plane of its 
rotation about the fulcrum /, to prevent vignetting and scattering of 
light from its walls. Suitable illumination is provided for the 
thermocouple junctions and for the scale of the air-mass meter. The 
thermocouple is connected with a D’Arsonval galvanometer whose 
deflections are read visually. When the knife switch, W, is reversed 
(only one side of the switch is used), the motor operates the crank 
and the image of the moon moves from one junction of the thermo- 
couple to the other, thus making the deflection. The snap switch, 
X, controls the microscope cover-glass, moving it into or out of the 
beam for the separation of planetary heat from the reflected sunlight. 
The pneumatic check, M, prevents jar, which might destroy the 
thermocouple or cover-glass when the magnet, N, acts. 

The apparatus was designed and constructed nearly ten years 
ago for use with a coelostat; but the coelostat appeared to introduce 
false deflections which could not easily be allowed for, particularly 
in daylight sky or when sunlight fell upon it. In the autumn of 1933 
it was attached to the 6-inch telescope mounting in order that it 
might be pointed directly at the moon and thus avoid these dis- 
turbances. 

OBSERVATIONS 

The observations consisted of (1) a set of deflections on the moon, 
five free and five with the cover-glass in place: (2) a set of five deflec- 
tions on a test lamp with the cover-glass in the beam. The test 
lamp was placed on the 60-foot tower telescope at a distance of 
71.32 m. The tower itself, however, gave a considerable deflection, 
owing to its self-radiation, which was included in the free deflection 
on the test lamp. The deflection with the cover-glass was not affected 
by this long-wave radiation, and ordinarily only this deflection was 
used, since the principal purpose of the comparisons was to test 
the constancy of the sensitivity of the apparatus from night to 
night. 

On part of the nights only one or two sets of deflections were made, 
when the moon was near meridian; but on many nights several were 
made at various air-masses to determine the atmospheric extinction. 
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On such occasions the air-mass meter was especially helpful, since 
the observations could be made for uniform intervals of air-mass, 
which greatly improves the plots from which the extinction is de- 
termined. 
REDUCTION OF THE OBSERVATIONS 
It has been pointed out’ that when we consider the radiation from 
a planet, the following relations hold: 


F=B+R, (1) 
Cg=BTr+Rtr , (2) 
_ Fir—Cg 
i= (3) 
o  RT ‘9—BT 
Rabat _ C8 _ (4) 


tr—Tr tr 


where F is the free deflection (the deflection without the cover-glass 
in the beam); Cg the deflection with the cover-glass; B and R the 
parts of F due to planetary heat and reflected sunlight, respectively ; 
Tr and ¢r the transmissions of the cover-glass for the planetary heat 
and the reflected light. It has also been shown that for moonlight 
the value of ér is very nearly o.go. 

The value of Tr as a function of the temperature of the planet 
has been tabulated,® but may be computed from the empirical 
formula 

Tr=4X1075(T— 264°)3/?, 264°<T< 700°. (5) 


For T < 264°, Tr=o. Below 700° K this equation gives the tabular 
value of Tr to the second decimal, which is sufficient for most pur- 
poses. We may estimate the maximum value of Tr, which occurs at 
full moon, from the following argument. The distribution of radia- 
tion over the full moon' is represented by 


E=a cos?/*6 , 


where @ is the angular distance of a point on the moon from the sub- 
solar point. Integrating this equation over the disk, we find that the 


6 Pettit and Nicholson, Mt. W. Contr., No. 392; Ap. J., 71, 102, Table I, 1930. 
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average energy is three-fourths the maximum value at the subsolar 
point. Since the temperature at the subsolar point is 374° K, the 
temperature corresponding to the average over the surface, by the 
fourth-power law, is 348° K, if we may regard the mixture of black- 
body energy emission as being itself a black-body emission. From 
equation (5) we see that 7r therefore has a maximum value of 0.031 
and is zero when T < 264°. Asa first approximation we may consider 
Tr equal to zero, whereupon equations (3) and (4) become 


Bo=F—1.11 Cg, (6) 
Ro=1.11 Cg. (7) 


The observations were all reduced with equations (6) and (7) and 
the mean curves corrected whenever 77 became a sensible quantity. 

It might be thought at first sight that equations (1) and (3) could 
be used to obtain R explicitly. This would necessitate a complete re- 
duction and discussion of all the data for B before the work could 
proceed, which would not be economical of the computer’s time; 
and, since the values of R, have already been formed in the pre- 
liminary determination of B, and B, their discussion can be carried 
on parallel with that of B,. The same explanation applies to the 
combination of equations (1) and (4). It was decided therefore to 
treat the two variables alike and to correct the mean curves for Tr 
in the final solution. 


ATMOSPHERIC EXTINCTION 

On a number of nights the reflected light and planetary heat were 
observed with the moon at various zenith distances. Plots of log R, 
and log B, against sec Z (air-mass) are shown in Figure 2 for Janu- 
ary 30. Although some curvature is noticeable, especially for the 
planetary heat, the plots are essentially straight lines. The loga- 
rithm of the transmission or extinction for unit air-mass can be ob- 
tained from these plots. For twenty-two nights on which suitable 
runs were made, the average transmissions for R, and B, were 0.90 
and 0.87, respectively, and the equivalent extinctions 0.11 and 0.15 
mag., respectively. 

A plot of the individual determinations against phase angle, while 








24 EDISON PETTIT 


showing some fairly definite connection, is discordant because posi- 
tive and negative phase angles do not agree in the general run, as 
might be expected since the equivalent temperatures of similar 
phases are substantially the same. The agreement is much better 
for reflected light than for planetary heat, however. 

The atmospheric absorption and the extinction as used here are 
distinct entities. The atmospheric absorption of planetary heat is 


r 


LOG DEFLECTION 
x 








7 5 
SECANT Z 


Fic. 2.—Plots to determine the atmospheric extinction. Planetary heat By and re- 
flected light R, separated by the cover-glass. 


nearly complete except in the: 8-14 » transmission band, while the 
effect of change in air-mass which measures extinction is principally 
confined to the wings of this band and is therefore much reduced. 
While the extinction amounts to only about 0.15 mag., the at- 
mospheric absorption exceeds 1 mag. 

Since in the case of reflected light from the moon we have a com- 
bination of scattering and absorption phenomena, the absorption 
and extinction are more nearly comparable, the absorption being 
about 0.4 mag. and the extinction 0.11 mag. 

On account of this relationship between absorption and extinction, 
the observations plotted in Figure 3 are those extrapolated to the 
zenith (unit air-mass) either by the observed extinctions for the night 
or by the mean value if only one or two observations were made. 
The upper curve (the heavy full line) in Figure 5 has been ob- 
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tained from the plot of B in Figure 3 by applying the absorption co- 
efficients of the atmosphere. 
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Fic. 3.—The full-line curves represent the reflected light Ro and the planetary heat 
By separated by the cover-glass of the thermoelectric lunar radiometer and are for the 
zenith on Mount Wilson. The dotted lines are Ro and Bo, corrected for the lunar plane- 
tary heat transmitted by the cover-glass. 


GEOMETRICAL CORRECTIONS 
The data from the ephemeris give the position of the moon rela- 
tive to the center of the earth. The value of sec Z for the position of 
the observer may be found by the formula 


: sin Z, 
sec Z=sec (Z+sin~ eee ) ; (8) 
V 3601—120 cos Zy 


where the last term is the parallax for the mean distance of the 
moon (60 terrestrial radii). For the purpose of computing the air- 
mass the empirical formulae 


log (A sec Z.)=1.7 log (sec Z>—1)—1.37 , sec Z,<8, (9) 
sec Z=sec Zo+A sec Zo, 
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are more convenient. The value of A sec Z, is 0.1 air-mass for sec Z, 
= 2.6, 0.2 for sec Z,= 3.5, and becomes 1.0 air-mass at sec Z,=7.3. 

The correction to sec Z, given by equation (g) applies to the center 
of the moon regarded as a disk. An additional, but much smaller, 
correction should be applied when the phase angle is large, since the 
actual center of gravity of the emitted lunar radiation is off center. 
If both rising and setting moon are observed to obtain the at- 
mospheric extinction, this correction will be nullified; and in the 
present work its value does not exceed 0°2, which might make a 
difference of 0.19 in the air-mass at sec Z,=7.3, where A sec Z, be- 
comes 1.0 air-mass. This correction has been neglected. Had sec Z, 
been linear in equation (9), it could have been neglected, as only the 
zenith values of R, and B, are used. 

Since the entire image of the moon falls on the receiver, a correc- 
tion must be made for the size of the image, which varies with 
diurnal parallax. To reduce the observed deflection to what it would 
be if observed at the center of the earth, we apply the correcting 
factor 

I 


~ 2 (10) 
30 sec Zo’ 


C=1.0003— 


which assumes the moon to be at a mean distance of 60 terrestrial 
radii. C varies between 1.0003 when the moon is on the rational 
horizon and 0.967 when in the zenith. 

A third correction must be applied to reduce all observations to 
mean lunar and solar distances. This factor was applied in the 
customary manner to the final values of B, and R, reduced to the 
zenith. 

PHASE ANGLE 


The phase angle is given by 
cos i=cos [180°—(M—S—E)] cos B , 


where 8 is the latitude of the moon, £ its elongation from the earth 
as seen from the sun, and M and S the longitudes of the moon and 
sun, respectively. The effect of 8 can usually be neglected except 
near full phase, and E cannot exceed g’. Since the plots of B, and R, 
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are made on co-ordinate paper divided into 2° units, the simpler 
formula 
i=180°—(M-—S) (11) 


was used for most cases. The error in (11) does not exceed 1° when 
+13°<i<—13°. 

During a series of measurements for atmospheric extinction, 
which may cover several hours, the phase angle changes. Since the 
plots are extrapolated to the zenith and the phase angle is computed 
for the moon at meridian transit, this change in 7 has no effect on 
the result. 

PLANETARY HEAT AND REFLECTED LIGHT 

The values of B, and R, plotted in Figure 3, which were obtained 
from the galvanometer deflections for each night, have been reduced 
to a standard test-lamp deflection of 59.7 mm. They are zenith 
values, reduced to the center of the earth, and correspond to the 
mean distances of the sun and moon. The heavy full lines represent 
the average values of these plots. We must now correct these curves 
for the planetary heat transmitted by the cover-glass, and in so 
doing obtain the curves of B and R. 

Substituting o.9 for /v in equations (3) and (4), we find the follow- 
ing simple relations 








Bo 
” gg, See isa) 
R=R,—1.11 BTr. (13) 


Since the temperature of the moon corresponding to average radia- 
tion from its surface is 348° at phase angle o°, we may obtain a first 
approximation for 7r by assuming that the plot of B, in Figure 3 
closely represents B and, applying the fourth-power law, obtain the 
temperature for each phase angle, and thus 77 from equation (5). 
It is found that Tr varies from 0.031 for i=o to zero for i= +72°. 
The broken lines in Figure 3 are the curves for B and R obtained in 
this manner from equations (12) and (13). They deviate so little 
from B, and R, that a second approximation is unnecessary. 

It will be noted that the corrected curves for both planetary heat 
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and reflected light are a maximum at full moon and sensibly sym- 
metrical about phase angle i=o°. At quarter phase the planetary 
heat has fallen to 17 per cent of the full-moon value, while the re- 
flected light is 15 per cent. The two corrected curves have been 
transferred to Figure 4 as heavy continuous lines. The light-line 
curves in this figure show the results obtained by Lord Rosse. The 
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Fic. 4.—Heavy full-line curves, R and B, observed at Mount Wilson; light full-line 
curves, observations by Lord Rosse. The broken line is R, observed by Stebbins and 
Brown with the selenium cell; the dotted line, R, by Rougier with the photo-electric cell 
and amplifier. 


broken line in the plot of reflected light is the curve obtained by 
Stebbins and Brown,’ and the dotted line, that of Rougier. The 
asymmetry found by some observers does not show to any great 
extent in the curves obtained with the thermocouple on Mount Wil- 
son, which in this respect agree with Rougier’s recent work. The 
ordinates of Lord Rosse’s curve of reflected light have been multi- 
plied by 1.11 to correct for reflection from the glass plate. Further 
corrections should be made for absorption in the glass, the thickness 
of which is unknown, for the reflecting power of the mirrors used 
(speculum), and for the greater atmospheric absorption at Birr 
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Castle. If these corrections were known, the curve would probably 
agree much better with that observed on Mount Wilson. 

Before discussing the appearance of these curves outside the at- 
mosphere, it will be convenient to replace the scale of galvanometer 
deflections by one of radiometric magnitude. It has already been 
shown! that the radiometric magnitude of planetary heat per square 
second of arc, m;+Am,;, is given by 


m;+Am;= 26.12—10 log T+Am,—2.5 log (1—1.11 Tr), (14) 


where Am; is a correction which makes the temperature of a point 
on the moon determined with this equation equal to that obtained 
from the amount of the absorbed solar radiation which can be ac- 
counted for. It is supposed that the existence of this correction is 
due largely to inaccurately known atmospheric absorption Am,. 
The values of m; and Am, as a function of T have been given in 
tabular form, but they may be represented by the simple equation 


m,+ Am; = — 2.6, 140° <T< 500° (15) 
within 5° of the tabular values of 7. Between 400° and 800° K the 
formula 

, ’ 1864 ) ) 
m,+ Am; = TP 73:35 400° < T < 800 (16) 


represents the higher temperatures very well, there being no sensible 
error up to 700° K. Formula (14), or the empirical formulae (15) and 
(16), give the radiometric magnitude of the planetary heat per 
square second of arc at the zenith on Mount Wilson as a function of 
absolute temperature. Further, the radiometric magnitude cor- 
responding to a 1 mm deflection may be determined by the use of a 
comparison standard. The radiation from a condensed filament 
Peerless Mazda lamp emitted at right angles to the filament was 
compared with that from a Hefner lamp in the laboratory at a dis- 
tance of 2 m, with the result 


g5 watt Peerless Mazda _ 


Hefner lamp 
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Since the Hefner lamp radiates 156X105 cal cm™ min™ at 1 m 
distance, the Mazda lamp radiates 9.02 10 * cal cm”? min, which 
is 83 per cent of the power input. Such a lamp working at 95 watts 
was set up beside the usual test lamp on the 60-foot tower telescope. 
After deducting the false deflection due to the tower, the value of A 
in the equation F = ACg was found to be 1.13. A test of the g5-watt 
lamp in the laboratory at a distance of only 56 cm, where no false 
deflection was observable, gave A =1.12. It thus appears that the 
procedure of subtracting the false deflection on the tower to obtain 
the free deflection on the light is correct, and that for the lamp at 
normal voltage we may obtain the free deflection from the equation 


F=1.12Cg. 


It was also found that the deflection for the free radiation from the 
95-watt lamp was equal to the radiation transmitted by the cover- 
glass from the test lamp multiplied by 1.165. Since all the observa- 
tions in Figure 3 have been reduced to a common test-lamp sensitiv- 
ity of 59.7 mm cover-glass deflection, the equivalent 95-watt lamp 
deflection for a distance of 7132 cm is 69.6 mm. The value £ in 
energy units for a 1-mm deflection is 

Bin RO ns eae ag calcm™ min“, (17) 

Rd 

where R is the distance of the lamp in meters and d the deflection cor- 
responding to 95 watts. The relationship between energy and radio- 
metric magnitude is 


E’!=17.3X10-"X2.512~"" cal cm™? min™ ; (18) 


hence the radiometric magnitude of lunar radiation for any deflec- 
tion in Figure 3 is given by 


m,= —26.905—2.5 log E’. (19) 


The scale of radiometric magnitude for the whole moon used for 
Figure 3 was obtained from equations (18) and (19). 
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The propriety of using the inverse-square law in equation (17), 
when the test lamp is at so great a distance, may be questioned. 
Calculation shows that for a distance of 71.32 m the infra-red ab- 
sorption bands of CO, would introduce a possible correction of 10 per 
cent. 

To test this matter, deflections on a lamp at this distance were 
compared with deflections on the same lamp at 2 m with the light 
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Fic. 5.—Planetary heat from the moon. Below, at the zenith on Mount Wilson; 
above, the full line is the observed curve for no atmosphere; the broken line, the theo- 
retical curve. 


cut down by a sector to 1 per cent. The sector was 16 inches in 
diameter and the opening was ruled on a circle dividing engine, then 
measured on the same instrument after cutting. In the zone em- 
ployed, the actual sector factor was found to be 0.0098. On one oc- 
casion a larger sector of factor 0.0500 was used. Tests made on three 
clear nights in Pasadena showed nearly io per cent on one occasion, 
5 per cent on another, and a null result on the third. The value 
obtained on Mount Wilson was 3 per cent. Altogether, it may be 
questioned whether this correction is real, and, if so, whether it is 
variable. It has not been used in the following computations. 

The lunar energy outside the atmosphere as shown in Figure 5 
may be obtained by subtracting Am, from the radiometric magni- 
tude of lunar planetary heat per square second of arc (m;+Am?}). 
This may be obtained from Figures 3 and 4 by adding the area of 
the moon in square seconds expressed in magnitude, 16.09. The 
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relation between these variables has already been given,° but, for 
convenience, is reproduced in Table I. 

The effect of Am,, the correction for atmospheric absorption, is to 
increase considerably the energy at crescent phases compared with 
that at full phase. 

It can easily be shown that if the moon were a smooth non-con- 
ducting sphere acting like a black body, receiving and emitting 


TABLE I 





m)+Am! Am, m.+Am! Am, 
TES SE eA eae aes | 6.35 mag. AMAR. 5 sas 1.74 mag. 
TORR eet ae chanics 5.08 eee 1.45 
SE Rene ar ee eee inary | acer 7 ie £230 
ten chstesaieea tates wienaacae Gs 3.20 Bi. 1.28 
ORO Rar eer arene nee ara 2.606 a I. 30 
5 2.16 SI hy pisces ite er 5.25 


energy proportional to the cross-section of the beam for every ele- 
ment of its surface, the phase-radiation relation for planetary heat 
would be given by 

— 4 Ry (r—1) cos i+sin 7 . fas 

3(1-+ 0s 12) T 

where &,, is the mean intensity of the planetary heat over the ap- 
parent surface projected on the celestial sphere, RX, the intensity at 
the subsolar point, and 7 the phase angle. Tables of the second fac- 
tor of this function have been given by Miiller.? The locus of points 
computed by this equation shown in Figure 5 is directly comparable 
with the observed curve for no atmosphere. 

The ratio of the planetary heat observed outside the atmosphere 
to that computed from equation (20) is unity when 7 =o, falls to 0.5 
when i= +75°, to 0.2 when i= +160’, and probably decreases up to 
i=180°. This large deviation of the observed from the computed 
intensity of planetary heat is possibly due, in part at least, to the 
occulting effect of the mountains on the very rough surface of the 
moon and, in part, to its transparency. 


7 Photometrie der Gestirne, Anhang, Table I, 1897. 
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RADIOMETRIC MAGNITUDE OF MOONLIGHT AND 
PLANETARY HEAT 

The radiometric magnitude of moonlight and of planetary heat 
may be read directly from Figure 3, or may be computed from the 
deflections. For full moon the deflections on reflected light and 
planetary heat, after correcting R, and B, for transmitted planetary 
heat by equations (12) and (13), are 154 and 584 mm, respectively. 
Multiplying by £’ in equation (17), we obtain for the intensity of 
radiation 


ER=39.2X107 cal cm™ min“, 


Ez=149X1077 cal cm™? min™. 


Applying equation (19), we find the radiometric magnitude of 
moonlight at full moon, M,z, and of planetary heat, M,,, to be 


M,r=—13.4, M,p=—14.8. 


The sum of Ep and E% gives for the radiometric magnitude of the 
whole radiation for the full moon 


Mpo=—15.0. 


In a previous paper’ it was shown that integrating the drift- 
curves across the full moon observed at the 1oo-inch telescope gave 


the value M,x=—13.3 mag. Since the corresponding mean phase 
angle was 12°, a correction of —o.19 mag. (Fig. 3) should be applied, 
making M,r;=—13.49 mag. The radiometric magnitude of the 


whole radiation of the full moon, according to these earlier data, 
would therefore be — 14.9. Considering the approximate nature of 
the integration of the curves and that the data applied to a specific 
case only, we may regard the agreement as fair. 

The visual magnitude of the full moon is usually® taken to be 
—12.55. If the heat index of moonlight can be determined, a value 
of the visual magnitude can be found from the present measures. 


§ Russell, Dugan, and Stewart, Astronomy, Vol. 1, Table IV, 1926. 
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The relationship of heat index to color index or to water-cell absorp- 
tion may be used for this purpose. 

It has been shown in a previous paper? that the color index of the 
moon on the Mount Wilson scale is 0.92 mag., corresponding to 
spectral class dKo, for which the observed heat index is 0.55 mag. 
The color index obtained from E. S. King’s photographic magnitude 
and the above visual magnitude is 1.18 mag., corresponding to dK6, 
approximately. This result gives 7/ =1.2 mag., a value almost cer- 
tainly too high. 

The water-cell absorption of moonlight may be obtained from the 
drift-curves' of lunar radiation made with the 1oo-inch telescope. 
If we correct the curve made with the cover-glass by equation (13), 
multiplied by 1.11, we find for the curves shown"? that WC =o0.61 
mag., which corresponds to H/ =o0.73 mag. in the table of observed 
coefficients for dwarf stars.” 

On February 22, 1923, a mountain just off the terminator on the 
dark side of the moon whose crest was illuminated by the sun was 
observed at the roo-inch telescope with the thermocouple” having a 
microscope cover-glass window. Since this mountain must have been 
at a very low temperature, no sensible part of the deflections can be 
attributed to planetary heat. The observation gave WC =0.63 mag., 
corresponding to HJ =o0.77 mag. on the dwarf branch. 

Since the color index is not a sensitive means of finding the heat 
index, the values obtained with it should perhaps receive but 
nominal weight. The true heat index is probably in the neighbor- 
hood of 0.75 mag. For this value we have for the visual magnitude 
of the moon 


M,=HI+M,r=— 12.63, 


in fair agreement with —12.55, the value usually given. Since there 
is at present little reason to suppose that H/ is a function of phase 
angle, the curve for M, should have the same form as that for M,,. 


9 Pub. A.S.P., 38, 242, 1920. 
Pettit and Nicholson, Mt. W. Contr., No. 392; Ap. J., 71, 102, Fig. 3, 1930. 
" Mt. W. Contr., No. 369; Ap. J., 68, 279, Table V, 1928. 


2S. B. Nicholson assisting. 
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TRANSMISSION OF LOW TEMPERATURE RADIATION 
THROUGH OZONE 

An absorption band of ozone occurs at 10 w in the middle of the 
8-14 w transmission band of water vapor through which planetary 
heat is measured. The thickness of atmospheric ozone is approxi- 
mately o.3 cm at N.T.P. and varies some 50 per cent. One may ask 
what effect this variation may have on the measured planetary heat 
from the moon. 

Since the temperature of the moon is in the neighborhood of that 
of boiling water, and since the radiation of the moon is measured 
with reference to space, we may simulate the conditions of observa- 
tion by measuring in the laboratory the transmission of ozone for 
radiation from boiling water relative to that of liquid air. 

To do this, a tube 130 cm long was closed at each end with rock- 
salt windows of 5 cm clear aperture and 6 mm thick. Inlet and exit 
stopcocks were attached near the ends. Ozone was produced in a 
Brodie (Berthelot) ozonizer fed with commercial oxygen, which was 
passed through a sulphuric acid dryer and a calcium chloride tower. 
A 50,000-volt transformer produced the silent discharge across the 
space in the ozonizer. Samples taken from the exit stopcock of the 
ozone tube and analyzed by the usual iodine method showed an 
ozone content of 2 cm at N.T.P. 

A blackened can of water, kept boiling by an electric heater im- 
mersed in it, was placed on a wooden slide beside a beaker of liquid 
air, a wood block separating the two. By moving the slide to and fro, 
the ozone tube could be exposed alternately to radiation from boiling 
water and from liquid air. The radiation passing through the ozone 
tube was received by a rock-salt lens of 1 inch aperture and 2 inches 
focal length, which formed an image of the blackened can or the 
beaker on a receiver of a thermocouple. 

Deflections were observed first with oxygen in the tube and then 
with ozone, after analysis. It was found that 93 per cent of the radia- 
tion passed through. While the ozone in the atmosphere is only one- 
sixth that used here, it is effective on the center of the 8-14 u trans- 
mission band of water vapor, which includes only one-fourth the 
radiation emitted by the moon. The absorption of lunar planetary 
heat due to atmospheric ozone is therefore about 4 per cent and may 





36 EDISON PETTIT 


vary from about 3 to 6 per cent. Probably the whole effect on a con- 
siderable number of lunar radiation measurements would not be 


serious. 


I am indebted to Dr. Fred E. Wright for measurements of lunar 
radiation throughout an entire cycle of phases. I am also indebted to 
Mr. Earl Karr and Mr. Glenn Moore for many of the measurements 
included in the program, and to Miss L. Ware for the reduction of 
the entire series. 

CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 
July 1934 





REMARKS ON LUNAR RAY CRATERS 
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ABSTRACT 

A short bibliography and résumé of various suggestions concerning the formation of 
lunar rays are given. The suggestion that the rays were formed by ejection from the as- 
sociated craters is examined quantitatively, and some aspects of the possible mode of 
formation are discussed. It is found that the minimum initial velocity required to pro- 
duce the longest ray from Tycho is 1 km/sec., taking into account the moon’s curvature. 
A “typical ray crater” has a (1) bright center floor, (2) dark annulus, (3) bright nimbus, 
(4) ray system. Assuming the projection to have taken place in the most efficient man- 
ner, the energy required to form the various features of Tycho is calculated (Table I). 
These energies are found to furnish no criterion of distinction between the volcanic and 
the meteoritic suggestions of formation of these craters, but are probably compatible 
with both. 


An equal-area map of the moon shows the distribution of ray- and light-surrounded 
craters listed by Elger. The distribution is analyzed by octants. There is a decided 
tendency for these craters to appear on the maria, probably owing to their better visi- 
bility there by reason of the contrast with the dark background. There is, however, no 
tendency toward grouping on the maria as between one octant and another. Similarly, 
the distribution off the maria is also random. The randomness of distribution favors the 
suggestion that the ray craters are of meteoritic origin. 

Any successful theory of the origin of the lunar rays must account 
for the following characteristics: (1) these rays are always associated 
with craters, and (2) emanate more or less radially from them; (3) 
they pass continuously over all sorts of country; (4) their visibility 
decreases with the angle between the line of sight from the earth and 
the direction of illumination from the sun (phase effect). (5) With 
the roo-inch telescope the rays are seen’ as the illuminated sides of 
low mounds. 

Many explanations of the rays have been attempted. Tomkins’ 
suggested that they might be alkali tracts similar to those which he 
found in India: one of these was goo miles long, 20 miles wide at the 
lower end, and 60 miles across at the other. He suggested that, owing 
to evaporation of moisture at some early time, alkali salts came to 
the lunar surface as white efflorescences taking the configuration of 
the country above the source of the material. In consequence of sub- 
sequent alternations of temperature the original alkaline crystals, in 
the absence of atmosphere, broke up into fine powder and remained 
visible. 

™ Pease, Pub. A.S.P., 36, 346, 1924. 2J.B.A.A., 18, 107 and 361, 1908. 
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Jeffreys’ at one time suggested that the rays are magma-filled 
cracks. The cracks were produced by tension due to cooling in the 
top layers of the lunar crust immediately following solidification. 
They grew downward and finally, by reason of the decrease in pres- 
sure, fusion of the lower rocks—-previously below their melting-point 

took place, and the magma was forced up the cracks by hydro- 
static pressure. This is practically the same suggestion that Na- 
smyth# also made. 

One might say that these are typical explanations of the origin of 
rays by surface and subsurface phenomena. Then there are the sug- 
gestions that the rays were formed by material falling down upon 
the surface. Chacornac seems to have been the first to suggest that 
the rays might be due to the deposition of fine white powder thrown 
out from the craters by volcanic action, upon a background of lower 
albedo. W. H. Pickering’ in 1892 recalled this explanation to mind. 
In 1919 airplane photographs of craters formed by exploding bombs 
dropped from planes led Ives® to believe that the lunar rays represent 
material ejected by meteoritic impact rather than by the volcanic 
action which Chacornac supposed. 

Buell and Stewart’ found that rays formed in the manner sug- 
gested by such an ejection theory would yield the phase effect. Their 
experimental study led to the conclusion that a lunar ray consists of 
a streak of scattered rocks of dark material matching in albedo the 
average surface of the moon, with the lighter-colored powder result- 
ing from the fragmentation of the rocks spread among these. It is 
assumed that this material was ejected from the associated craters 
either by volcanic action (Chacornac) or by meteoritic impact (Ives). 

It has recently been objected® that patches, rather than rays, 
would be formed by ejection of the material in this way—that the 
odds would be high against the formation of a ray by ‘‘one eruption 


3 The Earth (1924 ed), p. 147. 
1 Nasmyth and Carpenter, The Moon (London, 1874), p. 133. 
5 M.N., 53, 279, 1893. 


6 Ap. J., 50, 245, 1919. Attention is explicitly called to the fact that the photographs 
reproduced in this paper show clearly the actual formation of ray systems. One of the 
craters shown has a pronounced central peak. 


7 Pop. Astr., 40, 204, 1932. 8 J.B.A.A., 42, 308, 1932. 
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of a number of rocks, the individuals of which fall with remarkable 
uniformity of spacing and direction along a straight line, or a series 
of eruptions which no less remarkably distribute the rocks in this 
uniform way.” But at present we know little or nothing about 
splashes, and the only answer that can be given now is that the 
photographs by Ives (loc. cit.) show unequivocably the formation of 
rays.” 

It may be that at the time of the formation of a lunar ray crater 
the surface cracked radially in the immediate region around the 
crater and that these cracks served to guide the ejected material 
more or less uniformly within corresponding lateral limits. Until 
more becomes known about the configurations resulting from 
splashes, such criticisms and such answers to them as the foregoing 
can only be counted as personal opinion. 

Meanwhile, it is of some interest to examine quantitatively the 
suggestion that the rays were formed by ejection from the associated 
craters, and to study the distribution of ray craters on the lunar 
hemisphere. 

Taking into account the curvature of the moon’s surface, we first 
compute the minimum initial velocity at the most favorable angle 
of projection to form the longest ray, namely, the one from Tycho 
down across the Mare Serenitatis which covers a range of 78° (order 
of 1300 miles). This minimum velocity turns out to be only about 
1 km./sec. 

This value for the velocity does not exclude either the volcanic 
or the meteoritic hypothesis, being compatible with both. According 
to the meteoritic explanation, very much greater energies are avail- 
able, even exceeding the velocity of escape of 2.4 km/sec.’ In the 
case of terrestrial volcanoes" a velocity of 1 km/sec. is not impos- 

9 Whether or not “duds” with sufficient energy could produce the same effect as the 
explosive bombs is not known, and no attempt is made here to discuss this or to criticize 
Ives’s argument of the explosive action of lunar meteorites. See also F. R. Moulton, 
‘Some Quantitative Aspects of the Fall of Meteors,” Pub. A.A.S., 6, 372, 1931. 

© This, together with the fact that there are so few very long rays on the moon, 
might be submitted as a possible argument against the meteoritic theory. But there are 
too many unknown factors for definiteness pro or con. 

" Terrestrial volcanoes seem to be of two types: explosive, as at Lassen Peak, and 
lava-producing, as at Kilauea. The character of a given volcano, Day believes, depends 
on the gas content of the crystallizing magma. 
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sible. A block of rock was thrown from the Krakatoa crater’ for a 
distance of 30 miles and is estimated to have risen 31 miles in its 
flight. Neglecting air resistance and the size of the rock, v= 1 km/sec. 
and a=763°. In the case of Lassen Peak, California (May, 1915)," 
blocks of rock the size of a man’s hand were thrown for 1o miles, and 
horizontal blasts through the northeast side of the inclosing cone 
hurtled blocks of rock for more than 4 miles. Again v is of the same 
order. In general the maximum macroscopic velocity that a gas at a 
certain pressure and temperature can impart to a body ejected by it 
is of the order of the velocity of sound in the gas. If the gas is sub- 
ject to the adiabatic law,'+ then 


? 


ee eee ea (2-C r 
7 ~ Yo 4 Co— Cc’) 


where c is the velocity of sound and g the fluid velocity at pressure /, 
and where c, and gq, are the corresponding velocities at po. If p., which 
may be taken as the pressure in the discharging vessel, is very large, 
it may be legitimate to put g,=0; and as p diminishes, g increases, 
but only up to a certain maximum attained when p=o. If y=1.408, 
this limiting velocity is 2.214¢o. 

If Tycho, Olbers, and Anaxagoras may be said to represent the 
most general type, the “typical ray crater’’ may be defined as one 
having the following characteristics: (1) a bright crater floor; (2) a 
dark annulus outside and around the crater wall; (3) a bright nimbus 
beyond and around the annulus; (4) a ray system beginning at, and 
emanating from, this nimbus. Not all of the ray craters have all 
these features: Kepler, which is a notable example, has a bright floor 
but no dark annulus; the bright nimbus begins immediately outside 
the crater wall and stretches for 28 miles on the north and elsewhere 
from 50 to 70 miles, thereby covering about 9,000 square miles. The 
great ray system originates in this nimbus. Some craters have nimbi 
but no rays; these are called below “‘light-surrounded,” and are in- 


2 Krakatoa Committee Report, Royal Society, 1888; Iddings, Problems of Volcan- 
ism, Pp. 4. 
3 Day, J. Franklin Inst., 200, 161, 1925. 


4 Bull. Nat'l. Res. Council, “Report on Hydrodynamics,” No. 84, 529, 1931- 
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cluded in the class of ray craters. Still other craters have rays but 
no nimbi. 

On the present explanation the bright nimbi may be explained as 
being material thrown out in the same manner and of the same stuff 
as the rays but under a set of initial conditions more or less uniform. 
In some cases the prevailing initial conditions might have been such 
that there was a definite minimum range of projection of the same 
order as the width of the dark annulus. The bright crater floor can 
be explained as being due to material thrown back inside the crater; 


TABLE I 
TYCHO 
| | i 
| Assumed 
~~ Radius | Area } Volume Mean Initial | Energy 
‘“eature fe | . > . . > 
aaa (Mi.) | (Cm?) Cc) X10" | vx 104 | ErgsX 107? 
| Cm/Sec. | 
—_——- : 
2) TAPE We! 27 | 6X10%3 10.9 0.3 
: 
Wall Bh NAS Bed 0.5 | 0.01 
| | 
Annulus | 52 Ste Be) MeRapeN ret cry 
Nimbus go | 3X1035 2.4 5 6.0 
Rays | 3X10!'5 2.4 8 | 15.4 
| 
| | | 
. | | ~ > | 
Sum Re aden sia : Hake) Cae ee | 21.4 


this is a bit easier to visualize on the meteoritic explanation, the 
material being thrown inside from the edges of the impinging 
meteorite. 

Assuming the projection to have taken place in the most efficient 
manner, the energy required to form the various features of Tycho is 
given illustratively in Table I. The various radii are in miles; the 
other values are in c.g.s. units. Since the rays are seen with the 1oo- 
inch telescope at Mount Wilson as the illuminated sides of low 
mounds, their height is taken crudely as 25 feet. The density used is 
2, and the nimbus and the rays are taken to cover the same area and 
to be of the same volume. This may be too liberal an estimate for the 
rays, giving too large a value for the energy required to form them. 
The energy assigned to form the crater wall is no doubt much lower 
than was actually required: the value given represents merely the 
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potential energy of a wall of density 2, of height 6000 feet,'® and with 
both sides at a pitch of 60°. The great amount of energy expended 
in the rumpling of the wall and the surrounding region is neglected. 

Since no distinction is furnished between the meteoritic and vol- 
canic hypotheses by these energy relations, it is of interest to study 
the distribution of the ray craters. The distribution of lunar craters 
of all kinds is known not to be a random one. On the meteoritic ex- 








ic. 1.— Distribution of ray craters 


planation we certainly expect a random distribution. Figure 1 is an 
equal-area projection of the moon’s surface. The data are taken from 
Elger’s The Moon (1895) and are presented as an indicative, and not 
as an exhaustive, study of the distribution. 

The black dots represent light-surrounded craters; crosses alone 
are craters with rays but without nimbi; crossed dots are craters 
with nimbi and with rays. 

In Tables IL-IV the ray- and light-surrounded craters are not dis- 
tinguished from each other. Table II gives the distribution of these 
“ray craters” by octants with respect to the percentage of area of 

5 Schmidt, Charte der Gebirge des Mondes (Berlin, 1878), p. 290. The crater pit is 


also 6000 feet deep. 
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each octant covered by maria, as observed and as computed for a 
random distribution. 
The distribution in each octant with reference to the maria and 
non-maria is, therefore, not random: only one third of the lunar 
TABLE II 


DISTRIBUTION OF RAY CRATERS WITH REFERENCE TO MARIA 


Percentage of 


Area Obs. No. Obs. No. Comp. No. O-C 
Octant . 
Covered Or On On | On 
by Maria 
I 24 5 7 2.9 4.1 
2 47 7 3.8 3.2 
3 57 I S 2.3 0.7 
4. 71 fe) 10 hea 2.9 
5 38 3 ie O.1 0.9 
6 15 3 4 1.1 2.9 
7 fo) fo) ° 0.0 0.0 
8 21 4 8 2.5 5-5 
Sum 34 17 52 25.8 26.2 
TABLE III 
RANDOMNESS ON MARIA 
Percentage of 
Octant a se Obs. Comp. O-C 
— Total Maria — 
Bi Ash stoke aroun are ast 9 7 4.0 | 2.4 
2. 17 7 8.9 ce 
ae 21 3 10.9 | —7.9 
Aes 20 10 13.5 3.5 
oe 14 13 7 ieee 5.8 
6 5 4 2.8 22 
7 ° ° 0.0 0.0 
8 8 8 4.1 3.9 
: | 2 Bee 
SUM sos us Al 100 52 52.0 ro a Be 


hemisphere is covered by maria, but 52 of the 69 ray craters are on the 
maria. This apparent non-randomness of the distribution may be ex- 
plained: the bright rays and nimbi are more readily seen on the dark 
background of the maria than on the bright background of the non- 
maria. The seventh octant is noteworthy in this respect. Except for 
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a very small corner near the center, it is free of maria; and Elger lists 
no ray craters on it. The ray craters simply may be hard to dis- 
tinguish off the maria. Those that are seen there are either among 
the grandest systems on the moon or appear very faint. 

The question then arises whether the distribution of the ray cra- 
ters is random on the maria alone and random on the non-maria 


TABLE IV 
RANDOMNESS OFF MARIA 
Percentage of 


Octant om 5 Obs. Comp. O-C 
Total Non-Maria “_ 


i. 14 5 2A 2.6 
2. 10 I ay —O.7 
3 8 I 1.4 —=054 
4. 18) re) 0.9 —0O.9 
5. 12 3 2.0 [.O0 
a 10 3 2.9 0.3 
7 19 fe) ee. —3.2 
A ee 15 4 2.6 EA 

Sum. : . 100 17 10.9 O52 | 13 


alone. Tables III and IV give the answer as definitely affirmative. 
In so far as the number of ray craters listed in Elger may be taken as 
indicative of all ray craters that actually exist, the suggestion thus is 
favored that the ray craters are of meteoritic origin. 
PRINCETON UNIVERSITY OBSERVATORY 
PRINCETON, NEW JERSEY 
August 31, 1934 


THE SIGNIFICANCE OF RIGHT-ASCENSION ERRORS 
IN TRIGONOMETRIC PARALLAXES 
By T. E.STERNE 


ABSTRACT 


The significance of the data by van Maanen (A p. J., 78, 189, 1933), relating to sys- 
tematic errors, depending upon right ascension, in the trigonometric parallaxes of eight 
observatories, is examined by improved statistical methods. There is good evidence, 
from the harmonic analysis, for the existence of systematic differences, of period 24" in 
right ascension, between trigonometric parallaxes determined at the Allegheny and the 
McCormick observatories. The variation with right ascension cannot reasonably be at- 
tributed, entirely, to fluctuations of sampling. There is no evidence, from either the har- 
monic analysis or the tests for homogeneity, of any systematic variation with right as- 
cension, of the hourly mean differences between spectroscopic parallaxes and trigono- 
metric parallaxes determined at any of the following observatories: Mt. Wilson (60-in.), 
McCormick, Yerkes, Greenwich. A systematic error depending upon right ascension is 
suggested in the trigonometric parallaxes determined at the Allegheny Observatory. 
The errors, if present, are probably of period 24" in right ascension. Systematic errors 
of period 24" in right ascension are strongly suggested in the differences between the 
spectroscopic and the trigonometric parallaxes determined at the Sproul and the Yale 
(photographic) observatories, both by tests for homogeneity and by harmonic analysis. 
There is some evidence for the existence of a systematic variation, of period 6" or less in 
the right ascension (but not necessarily exactly periodic), in the differences between the 
Cape photographic and spectroscopic parallaxes (the variations of longer period can be 
attributed to fluctuations of sampling). The point is not stressed, despite a strong indi- 
cation of inhomogeneity, because of the inherent implausibility of such short-period 
variations. There is no evidence for the existence of systematic error depending upon 
right ascension, in the spectroscopic parallaxes used in the comparisons. Hence the 
right-ascension errors, in the few cases mentioned above where they are probably pres- 
ent, are most probably in the trigonometric parallaxes. 


Following an earlier paper,’ van Maanen has published a more 
complete compilation’ of data relating to what he claims are system- 
atic errors, depending upon right ascension, in the trigonometric 
parallaxes of stars studied at eight different observatories. In Table 
I of the latter paper he classifies the stars that have been measured 
for parallax both at the Allegheny and at the McCormick observa- 
tories into twenty-four classes, according to the hours of their right 
ascensions, and lists the number of stars, as well as the mean dif- 
ference + Allegheny — 7 McCormick, for each class. In his Figure 1 
he plots, against each hour of right ascension, the average of the 
hourly mean differences for the hour in question, the hour preced- 


t Mt. W. Contr., No. 357, 1928. 
2 Ap. J., 78, 189, 1933; Mt. W. Conir., No. 474, 1933. 
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ing, and the hour following. The figure suggests at first sight a 
systematic dependence of the three-hour means upon right ascen- 
sion, and van Maanen finds that the plot can be fitted fairly closely 
by a trigonometric series, in right ascension, involving seven con- 
stants which he obtained by the method of least squares. He con- 
cludes that there can be no doubt that these differences are a 
function of right ascension. He likewise treats, for eight observa- 
tories separately, the overlapping three-hour means representing the 
differences between the trigonometric parallaxes and a common sys- 
tem of spectroscopic parallaxes selected from the list of W. S. Adams. 
The plots are shown in van Maanen’s Figure 2, and the hourly means 
and numbers of stars in each hour group are given in his Table V. 
The plots all suggest, at first sight, a systematic dependence of the 
overlapping three-hour means upon right ascension, and each plot 
can be fitted fairly well by a trigonometric series involving seven 
constants, peculiar to the plot, obtained by least squares. Reasoning 
that the common system of spectroscopic parallaxes is not likely to 
contain errors depending upon right ascension, he concludes that 
there are systematic errors in all the trigonometric parallaxes, de- 
pending upon right ascension, and advocates their removal through 
the use of the corrections obtained from his least-squares solutions. 
Between successive, overlapping, three-hour means of the sort 
plotted and analyzed by van Maanen, there must be a correlation of 
+, even when the successive hourly means are not correlated. It is 
the correlation between successive points which suggests to the eye 
their systematic dependence upon right ascension; and after the high 
correlation of +3 is artificially introduced, even a set of quite ran- 
dom residuals will show the sort of irregularities apparent in van 
Maanen’s Figures 1 and 2. To this extent the use in his figures of 
three-hour overlapping means was unfortunate. I show in Figure 1 
the plots of the separate hourly means, drawn from van Maanen’s 
Tables I and V; they show no very striking correlation with right 
ascension. The irregular dispersions are caused largely by the un- 
equal weights of the points, since in some hour classes there are more 
stars than in others. The method of least squares automatically se- 
lects the values of the seven constants that give the very best fit 
possible between the trigonometric series and the set of three-hour 
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overlapping means. With seven disposable constants, it is entirely 
possible that the smoothed data? shown in van Maanen’s Figures 1 
and 2 could be fitted as well as they have been fitted or better, even 
if the original hourly means were not correlated at all with right 
ascension, but merely fluctuated about a mean value in an entirely 
random fashion. 
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I'1G. 1.—Mean differences 7,—7, grouped according to hours of right ascension 


Since it is important to answer the question as to whether or not 
there are systematic errors in the trigonometric parallaxes, depending 
upon right ascension (particularly in view of van Maanen’s sugges- 
tion that corrections, depending upon right ascension, should be ap- 
plied? to the parallaxes determined at the various observatories), we 


3 Only eight of the points in each of his figures are independent. 


4In Table IX he shows that the use of his system of corrections does on the whole 
reduce the discrepancy between the trigonometric parallaxes determined at different 
observatories; but the same reduction would have been effected whether the errors, for 
which he applied corrections, had been systematic or random. To correct for random 
errors by systematic corrections depending upon right ascension would of course not be 
legitimate. It is only after errors are known with certainty to be systematic that sys- 
tematic corrections should be introduced. 
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present what is intended and believed to be a rigorous statistical dis- 
cussion of the material in van Maanen’s Tables I and V. We use two 
different methods of investigation, to which the material in the 
above-mentioned tables is excellently adapted.s 


HARMONIC ANALYSIS 


Van Maanen solved each of his sets of overlapping three-hour 
means by least squares, using twenty-four equations of condition, 
for a constant term and for sine and cosine terms having periods of 
twenty-four, twelve, and six hours in right ascension. He did not 
publish the mean errors of the resulting constants necessary for tests 
of significance. We have therefore repeated his analysis, with the 
difference that we have (1) used hourly means instead of overlapping 
three-hour means in the equations of condition, and that we have 
(2) sought to determine only three constants instead of seven, by 
confining our attention to the constant term and to terms having a 
period of twenty-four hours in right ascension—terms inherently 
more likely, perhaps, than those of shorter period. With each hourly 
mean used as the second member, M/, of an equation of the form 


a+b sin g+ccos g=M , 


the constants a, b, and c were determined by a least-squares solution. 
The weight of each equation was taken equal to the number of stars 
involved in its determination. The values of the resulting constants, 
with their mean errors, all in thousandths of a second of arc, are as 
shown in Table I. 

We test the significance of the results for b and c, leaving the dis- 
cussion of the significance of the a’s to those who are familiar with 
the possible sources of constant systematic error. It would be legiti- 
mate to test the significance of the 6’s and c’s separately, since the 
correlation between } and c found from the least-squares solution is 
in all cases small, ranging between +0.07 and —o.07 in the nine 

5 See also C. Davidson, ‘Systematic Errors in Trigonometric Stellar Parallaxes,” Ob- 
servatory, 57, 236, 1934, for further comments on van Maanen’s paper. Note added 
October 24: Fora rediscussion of the material following the methods of van Maanen, the 
reader is referred to a recent paper by S. A. Mitchellin Ap. J., 80, 200, 1934. I acknowl- 
edge gratefully the opportunity of reading Professor’s Mitchell’s paper in proof. 
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cases above.® From tables of the normal error function, the proba- 
bilities could be found of obtaining, through fluctuations of sampling, 
values of b and c numerically greater than those actually found. But 
we are more interested in the significance of b and c jointly. Neither 
the amplitude (?+c’)?, nor the amplitude squared is, even roughly, 
a normal deviate, although } and ¢ aré normal deviates. We ignore 


TABLE I 


Series | a b c 
aw Alleg.—a McC............ | —1.02+0.88 | 3-59t3.13 | ~—1.99¢1.33 
x Mt. W. 60—-2 Spec....... | 3.78+1.42 | 0.83+1.95 —2.64+2.04 
2: AlGE. ~ 8 SpOC., 0. -2: +40 ~COELO 78 | 2.21+1.03 | —1.03+1.27 
w MeC.—x Spee... .2.6525+}| =O20EO75 | —O.007F.00- | 0O.43+1.11 
wx Yerkes—a Spec........... 2: 2E 21s | 2.24+2.7% | —2.443.14 
a Sproul—-z Spec........ 4.88+1.56 | 7 ISEG.33 | =—0.4922.42 
w Greenw.—7 Spec......:...| =On7IED.53. | —P.88E2i14 | —320042.57 
aw Yale Phot.—7m Spec.......| 6.38+1.43 —4.45+1.92 | —5.37+2.05 
a Cape Phot.—7z Spec...... . | 3.46+2.22 0.56+3.17. | —6.25+3.12 


the trivial correlation between } and c. It can easily be shown that 
if t, is the ratio of b to its mean error, and /, that of c to its mean 
error (/, and ¢, being in fact almost normally distributed, with unit 
standard deviation), the probability that a value of ¢+é? will be 
obtained, through fluctuations of sampling, larger than or equal to 
the value observed, is e-($+)/2.. This probability, called P,,, is listed 
below; it tests the significance of b and c jointly. 


Series Lie Series Pee 
mr Alleg.— a McC....... 0.004 2#Sproul—zSpec....... 0.003 
x Mt.W. 60—7 Spec.... .39 mw Greenw.— 7 Spec......  .16 
w Alleg.— w Spec........ .07 x Yale Phot.—x Spec... .002 
w McC.— = Spec........ 62 a Cape Phot.— 7 Spec... 0.13 
wx Yerkes— x Spec....... 0.53 


The ordinary rule in astronomical work has been to regard a nor- 
mal deviate as significant when it exceeds three mean errors, for 
which the probability of obtaining that deviate or a numerically 
larger one of the same sign is but 0.13 per cent; 1 per cent, corre- 

6 The coefficient of correlation is in each case the quotient of the co-factor of [bc] in 
the determinant of the coefficients in the normal equations, divided by the geometrical 
mean of the co-factors of [bb] and [cc], and would vanish completely but for the unequal 
weights of the 7’s. 
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sponding to 2.33 mean errors, is more lenient and probably strict 
enough in practice; while 5 per cent, corresponding to 1.65 mean er- 
rors, is decidedly lax. The value of P for m Alleg.—aMcC. falls 
short of the 1 per cent point and must be considered significant; a 
real error depending upon right ascension is very probably present 
in one or both sets of trigonometric parallaxes. The sine-curve in 
the  Mt.W. 60—7 Spec. comparison cannot be said, with any as- 
surance, to be more than a mere statistical fluctuation; nor can the 
sine-curves for McCormick, Yerkes, Greenwich, and Cape Photo- 
graphic. Since the McC. —z Spec. results are not significant, while 
the a Alleg.—a McC. are significant, the cause is perhaps to be 
found in a right-ascension error in the Allegheny parallaxes, a sug- 
gestion that is considerably strengthened by the value 0.07, in itself 
not particularly significant, for the m Alleg.—z Spec. comparison. 
The sine-curves for Sproul and Yale Photographic (in Johannesburg) 
are very probably significant. The significant sine-curves are drawn 
solidly in Figure 1, the curve for Allegheny, of more doubtful sig- 
nificance, being broken; none of the others have been drawn. The 
values of P listed above, in the absence of systematic errors having a 
period of twenty-four hours, should be uniformly distributed, at ran- 
dom, between o and 1. After the values of P for m Alleg.—a Mc.C., 
a Alleg.— a Spec., 7 Sproul—7z Spec., and mw Yale Phot.—~7 Spec. 
are removed from the list, the remaining values cannot be said, col- 
lectively, to be significant. Two of the values are larger than 0.50, 
and three are smaller. The observed mean is 0.37, with an expected 
mean error of 0.13 equal exactly to the amount by which the ob- 
served mean falls short of the expected mean. 


TESTS FOR HOMOGENEITY 

In the preceding section we computed by least squares particular 
parameters, namely, the coefficients of the sine and cosine terms, 
and tested them for significance. There would have been no par- 
ticular reason for testing these parameters rather than any others 
of the enormous set of possible parameters, had it not been for the 
strong presumption that if any systematic errors were present de- 
pending upon right ascension, they would be seasonal, and with a 
period of twenty-four hours in right ascension. To test for the reality 








ERRORS IN TRIGONOMETRIC PARALLAXES 51 


of a sine-curve of a particular period is artificial; and testing for the 
reality of any sine-curve is somewhat artificial, for even if the sus- 
pected periodicity is actually present, it is likely that more than a 
single term in a Fourier expansion will be needed to represent it. 
There is fortunately a simple and elegant general method, devised 
by Professor R. A. Fisher,’ for dealing with such problems; he calls 
it the ‘‘analysis of variance.”’ It does not depend upon any arbitrary 


TABLE II 











—$_—_— j 7 

Hours | oh+ | 4b+ 8h+ 12h-+ | 165+ 20h-+ 

é = Ss oe | me 5 pee ten ae 
Oi. 254.13: | 623 52} 350-14) = 525 S| 1.0 25 | —8.3 33 
I —8.3 18 7.3 16 24 FA] KO oO | —4.7 28 | —8.8 22 
2. —2.3 18 O.2 IS: =O. E FE},  ‘O-5 45 | -— F-02595 0.5 25 
3 5-113 | 7-418 | —3.1 8 | —3.1 21 | —3.3 45 | —4.1 15 

Means..| —1.50 4.47 —o.0o1 | —3.60 | —2.60 —5.44 

Weights. 62 64 47 51 | 123 CO 95 

| 





assumptions as to the exact nature of the inhomogeneity of the ma- 
terial tested. An application of it to van Maanen’s data will supple- 
ment the investigations of the previous section, and will also pro- 
vide an illustration of its usefulness in problems arising in astronomy. 

We cast the data found in van Maanen’s Table I into the form 
shown in our Table II. The hourly means have been grouped into 
six classes of four hours each. Each entry consists of two parts. 
First there is given the mean difference of trigonometric parallax, 
a Alleg.—a McC., for the stars in that particular hour group that 
have been measured at both observatories. Next is given the num- 
ber of stars in the group, hereafter called the ‘“weight.”” Below each 
four-hour class is given the weighted mean difference for the class 
(hereafter called the “‘mean of the class’’) and the total weight for 
the class. The sum of all the weights is 442, and the weighted mean 
of all the classes is —1.87. The units used are thousandths of a sec- 
ond of arc. We square each hourly difference, multiply it by its 
weight, and sum over all the 24 hours. The result is 11,674. We 
square the mean of each class, multiply it by its weight, and sum 
over all the six classes. The result is 5721. Finally, we square the 


7 Statistical Methods for Research Workers (Oliver & Boyd, 1932), chap. vii. 
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On 
No 


weighted mean of all the classes and multiply by its weight, obtain- 
ing the result 1549. In Table IIa, column 3, the top entry is the re- 
mainder 11,674~5721; the middle entry is 5721-1549; the bottom 
entry is 11,674-1549. The reader will recognize these numbers as 
merely the weighted sum of the squares of the individual entries in 
Table II about the weighted means of classes, the weighted sum of 
squares of means of classes about the weighted mean of all the 
classes, and the weighted sum of the squares of the individual entries 


TABLE Ila 


aw ALLEG.—a MCC. 


ene? Degrees of | Weighted Sum | Mean : 
Variance : ; ; S.D. 
Freedom | of Squares Square 
Within classes............ 18 5953 331 | 18.2 
Between classes............| 5 4172 834 28.9 
jo eae eee ee 23 10125 | (s=0. 4064) 


4" classes. Inhomogeneity not suggested. 


about the weighted mean of all the classes, respectively. ‘lhe total 
number of ‘“‘degrees of freedom” within classes is 18, since if the 
mean of a class is fixed, only three entries can vary independently 
in the class. Since the mean of all the classes is fixed, there are but 5 
degrees of freedom between means of classes. The total, 23, is the 
number of degrees of freedom for all the hour groups if the mean of 
all the classes is fixed. Dividing the sums of the squares (in the third 
column) by the degrees of freedom, we obtain the entries in the 
fourth column; their square roots are given in the last column. It 
will be recognized that these are estimates of the standard deviation 
of an hourly mean of unit weight, the value 18.2 being an estimate 
obtained from within classes, and the value 28.9 being an estimate 
obtained from between classes. In the absence of systematic differ- 
ences between the six classes, these estimates should be equal, except 
for possible differences due to mere fluctuations of sampling in the 
data. A significant excess of the “external” estimate over the ‘‘in- 
ternal” estimate would prove that the classes differed significantly 
from each other and hence that there were systematic errors depend- 
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ing upon right ascension in the trigonometric parallaxes obtained at 
one, or both, observatories. 

It can be shown that if the classes are homogeneous, the two esti- 
mates of the standard deviation behave exactly as though they had 
been obtained from two samples, one of nineteen individuals and 
the other of six individuals, drawn independently and at random 
from the same universe, of unknown standard deviation o. The sta- 
tistical distribution of v, the ratio of the ‘external’ estimate of the 
standard deviation to the “internal,” is for a normal universe exactly 


Nita 
N2 


p(v)dv=const. Xv" (a+ ‘| : dv , 


n 


where 7, and n, are the external and internal degrees of freedom, 
respectively. When m, and m, are nearly equal, then the function 
p(v) is nearly symmetrical about the value 1, and z=log, 7 is dis- 
tributed about the mean value o in a nearly Gaussian fashion, with 
a standard deviation given roughly by [(7,+2.)/2n,n,]'. Fisher (loc. 
cit.) has given two tables of the exact values of z (corresponding to 
various values of m, and m,) such that the chances of exceeding them 
are 5 and 1 per cent, respectively. Of course, the 5 and 1 per cent 
points for negative values of s (i.e., values such that the chances that 
z will fall short of them algebraically are 5 and 1 per cent, respective- 
ly) are obtained by interchanging m, and , as the arguments of the 
tables. In our case, a large negative value of z would indicate the 
opposite of a dependence upon right ascension. From Table Ia, we 
find that s=0.464. For n,=5, n,=18, the positive 5 and 1 per cent 
points are 0.5099 and 0.7232, respectively, while the negative 5 and 
I per cent points are —o.7605 and — 1.1312, respectively. The value 
of s that we have obtained from the data of van Maanen’s Table I 
does not even reach the 5 per cent point, and cannot be regarded as 
significant. There is thus no evidence of any systematic dependence 
upon right ascension, of the differences between the Allegheny and 
the McCormick trigonometric parallaxes. Such slight inhomogeneity 
as there is could easily have arisen as a mere fluctuation of sampling. 

There are three considerations governing the size of the classes 
into which we should divide the hourly mean differences. First, un- 
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less the classes are small, there is a chance that such errors as may 
be present, of short period in right ascension, might be smoothed 
out and overlooked. Second, by dividing the hourly differences into 
twelve classes of two hours each, instead of into six classes of four 
hours each, the degrees of freedom would become 12 and 11 instead 
of 18 and 5; resulting in comparable accuracies for the internal and 
external estimates of the standard deviation, in a very nearly normal 
distribution of z, and in the possibility of easily assigning actual 
probabilities to the values of z found. These two considerations both 
point to the advantage of using two-hour, instead of larger, classes. 
On the other hand, there is a third consideration, which assigns an 
advantage to larger classes. The weights of the hourly differences 
are comparable; let us suppose for simplicity that they are all equal. 
If there is present a real error depending upon right ascension, con- 
sisting of periodic terms the sum of whose amplitudes squared is A’, 
and if we use classes of V hours each, then unless the classes are so 
large that some of the terms are smoothed out, the (standard de- 
viation)? of the means of classes will be nearly }A?+a?/N, where a is 
the standard deviation of the random errors that are superposed on 
the systematic ones. We shall thus be led to an external estimate of 
the (standard deviation)? nearly equal to }A*V +a’. The internal 
estimate of the (standard deviation)’ will be very nearly a’. The 
ratio v? will thus be nearly equal to 1+3A?N/a’, on which of course 
will be superposed the errors of sampling. If systematic errors are 
present, we shall therefore increase the chance of discovering them 
by using larger classes, provided always that we do not smooth out 
the periodic terms by so doing. It may be shown that for discovering 
periodic errors of a period of twenty-four hours, of an amplitude 
which is small in comparison with the random error a, classes of 
eight hours each should yield the largest ratio of the external to the 
internal estimate, smoothing being taken into account. For a period 
of twelve hours, four-hour classes should be used; and for a period of 
six hours, two-hour classes. The results obtained by dividing the 
data of van Maanen’s Table I into twelve classes of two hours each 
are as in Table IIb. The value of z is 0.358; the positive 5 and 1 
per cent points are 0.4994 and 0.7194, respectively; the negative, 
—o.5126 and —o.7405. The inhomogeneity is therefore not signifi- 
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cant. Finally, we repeat the analysis, using three classes of eight 
hours each. The results are as in Table IIc. The value of z is thus 
0.565. The positive 5 and 1 per cent points are 0.6216 and 0.8772; 
the negative, — 1.4839 and — 2.2999, respectively. Thus there is no 
good evidence for inhomogeneity in the data of van Maanen’s Table 
I, no matter how we treat it. Good evidence was obtained, however, 


TABLE II} 


wm ALLEG.—a MCC. 








l | : 
| | 
Sa | Degrees of | Weighted Sum Mean 2 
Variance | : | : | ; | S.D 
| Freedom | _ of Squares | Square 
SAS é aay] en eer mt ae 
Within classes......... 2 12 | 3524 204 | yen 
Between classes......... a 6601 | 600 24.5 
: : : “ail fee Say eee 
[SOE as ene ee 2 | 23 | 10125 livres xen | (s=0.358) 
2 classes. Inhomogeneity not suggested. 
TABLE IIc 
x ALLEG.—a MCC. 
Ge Degrees of | Weighted Sum | Mean . 
Variance | : | : S.D 
Freedom | of Squares | Square 
Within classes............ | 21 | 7820 | 7 | 19.3 
Between classes............ 2 | 2305 | 1152 | 33-9 
| ee ee 23 | 10125 | EER Oa aes | (s=0. 56s) 
| | 


| 





8 classes. Inhomogeneity not suggested. 


for the existence of a right-ascension error, from the harmonic analy- 
sis of the preceding section; and we must still consider right-ascen- 
sion error to be present, despite the negative results of the search 
for statistical inhomogeneity. 

The same procedure as that by which we analyzed van Maanen’s 
Table I has been applied to his Table V, containing the mean hourly 
differences, for each of eight observatories, between the trigonomet- 
ric and spectroscopic parallaxes. The tables of results, analogous to 
our Table Ila, follow (Table III). The values of z are tabulated in 
Table IV. 
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TABLE III 
a MT. W.—7z SPEC. 


eat Degrees of 
Variance 


Freedom | of Squares 
Within classes........ 12 | 3658 
Between classes............| II | 4767 
Total | 23 | 8425 


2" classes. Inhomogeneity not suggested. 


a ALLEG.—7 SPEC. 


ance Degrees of Weighted Sum 
Variance * a 
Freedom | of Squares 


Within classes ee 12 61609 
Between classes. ...... =") II 6247 
MOUS Coe tee chain siinacs 23 


| 
12410 


2" classes. Inhomogeneity not suggested. 


aw MCC.—-x SPEC. 


a Degrees of 
Variance 


Freedom of Squares 
Within classes. ..... bs 12 4478 
Between classes............] II 4194 
NAGA ec nce acise es 23 8672 


2" classes. Inhomogeneity not suggested. 


mw YERKES—z7 SPEC. 


Degrees of Weighted Sum 


Variance ; . 
Freedom of Squares 
Within classes... 52... | 12 11436 
Between Cineses:...........:] II 7927 
1 nS ee eee SES cee 23 19363 


25 classes. Inhomogeneity not suggested. 


| Weighted Sum | 


| Weighted Sum | 


Mean 
Square 
| aa 
| 305 
| 433 


Mean 
Square 


514 
568 


Mean 
| Square 


| 373 
381 


Mean 
Square 


S:D; 
E75 
20.8 


(s=0.176) | 


NO ND 
wn 
on 


(s=0.047) 


} 

S:D. 
| 

19.3 

19.5 
3=0.010) 
, 
| 

$.D. 
39.9 ‘ 
20.9 ) 

$= —0.139) 

j 
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TABLE IlIl—Continued 


7 SPROUL—7 SPEC. 











_ 
ae Degrees of Weighted Sum Mean - 
Variance : | _ | 3 $.D 
| Freedom | of Squares Square 
= MAES eee eaten & ere SS 
| | 
Within classes i... cos elcnes | 12 5910 492 22.2 
Between classes........... | II 11432 | 1039 32.2 
i a Ce See ST a 
oS ete ae 23 17342 | SNS nee (s=0.372) 
| | 
2" classes. Inhomogeneity not suggested. 
a GREENW.—7 SPEC. 
Be 
Widaeeds Degrees of | Weighted Sum Mean SD. 
} Freedom of Squares Square 
- | | - - | — —_ _ —_-—— —_—__— $$ 
Within classes........ Sahesd 12 5106 426 20.6 
Between classes............| 11 4172 379 19.5 
= | . ell I ee ne a ee 
Ee Lonny | 23 | 9278 reer se t= —0.055) 
} 2" classes. Inhomogeneity not suggested. 
x YALE PHOT.—7 SPEC. 
| > no 
P | Jegrees | Weighted S Mez 
ee on I egrees of eightec um | fean S.D. 
| | Freedom | of Squares | Square 
Within classés..........5.... 12 5535 461 21.5 
Between classes........... | II 10392 045 30.7 
— Cech seems ee TS 
Co eee 23 15927 | ee ee | (s=0.357) 
| | 
, k . 
2" classes. Inhomogeneity not suggested. 
a CAPE PHOT.—7z SPEC. 
| = = A TET 
Degrees Veighted S | N | 
ee egrees of \ eighted um fean SD. 
Freedom | of Squares Square | 
, Within classes...0..400..0. 12 2043 | 170 13.0 
Between classes............| Il 8212 | 747 27.4 
at — Co ia aia ii 
Bettie oo oc eee 23 10255 | center year | Z=0.747) 
| 





25 classes. Inhomogeneity strongly suggested. 
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TABLE I1I]—Continued 
aw MT.W.—7 SPEC. 


| } 
ss faee | Degrees of Weighted Sum Mean ; 
Variance : ‘ . S.B. 
Freedom | of Squares | Square 
| | 
ye . TP en) } > we wmQ | | re 
Within classes............. | 18 557 | 310 | 17.6 
Between classes........... 5 2847 569 23.8 
| 
- - | —_— —— | —— - 
= . | 
LGHAl. ee aees cae 23 | 8425 | er ee | (s=0. 300) 


4" classes. Inhomogeneity not suggested. 


a ALLEG.—7z SPEC. 


Degrees of | Weighted Sum 


| | 
| Me: 
ee a | swlean S.D 
Variance Freedom of Squares Square : 
Within classes....... el 18 9173 510 22.6 
Between classes........... 5 3243 | 649 | 25.4 
ea i | ; 
eC RE ae eae Mee i 23 | 12416 Jesse eeeee ees | (s=0.113) 
mig | | i | 
4" classes. Inhomogeneity not suggested. 
a MCC.—7z SPEC. 
| | 
: Degrees of | Weighted Sum | Mean ‘ 
Variance | . ; . S.D. 
Cees | Freedom | of Squares Square | 
| 
Within classes........ ciel 18 6245 | 347 18.7 
Between classes............| 5 | 2427 | 485 22.0 
MOG eee nc. a 23 | 8672 | epee ae ie | (s=0. 165) 
| 
4" classes. Inhomogeneity not suggested. 
a YERKES—z SPEC. 
| : | 
ee Degrees of Weighted Sum | Mean | : 
Variance : . | ; S.D. 
Freedom of Squares Square | 
MVACIANICIASSES 5 5.6 os 6 os 3s a 18 | 15522 862 | 20.4 
Between classes............| 5 | 3841 768 | a7 7 
2 aera gers | 23 a SO eee ...|(2=—0.058) 


| 


4" classes. Inhomogeneity not suggested. 
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Variance 
Within classes............ 
Between classes. ... 
Total 


| 


| 
| 
| 
| 


TABLE II1]—Continued 


ma SPROUL—7 SPEC. 


| Weighted Sum 
of Squares 


Degrees of 
Freedom 


18 9304 
5 | 7977 
23 17342 


4" classes. Inhomogeneity suggested. 


Variance 

Within classes 
Between classes. .... 
fis: | ee 


a GREENW.—z7 SPEC. 


| Weighted Sum 
of Squares 


Degrees of 


Freedom 


18 7025 
5 2253 
23 9278 


4" classes. Inhomogeneity not suggested. 


Variance 


Within classes 
Between classes........ 


Total: .....; 


a YALE PHOT.—7z SPEC. 


n 
| 


} 
| Weighted Sum 
of Squares 


Degrees of 
Freedom 


18 7809 
5 811d 
23 15927 


4" classes. Inhomogeneity suggested. 


Variance 
Within classes. . 
Between classes 


Total 


a CAPE PHOT.—7z SPEC. 


Weighted Sum 
of Squares 


Degrees of 
Freedom 


18 | 8076 
5 2179 
23 | 10255 


| 


4" classes. Inhomogeneity not suggested. 


Mean 
Square 


520 
1595 | 


| 
Mean | 
Square | 


300 | 
451 | 


Mean 
Square 

434 

1624 





Mean 
Square 


449 
436 
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SD. 


21.2 
20.9 


in —avnnal 
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TABLE II[—Continued 
a MT.W.—7z SPEC. 















































pie Degrees of | Weighted Sum | Mean 
Variance . . S,D. 
Freedom of Squares | Square 
| 
SERS ; Aa ; —_ ; ; | ’ 
Within classes............. 21 7513 | 358 18.9 
Between classes............ 2 | QI2 | 450 21.4 
1S RP ah nae a 23 | 8425 | ROSE EE a Th | (g=0.122) 
} 
8 classes. Inhomogeneity not suggested. 
a ALLEG.—7z SPEC. 
| | | 
ade Degrees of | Weighted Sum | Mean | . 
Variance : : | ; S$.D 
Freedom of Squares | Square 
Within classes: .....0:5 6.50004. 21 9528 454 | 21.3 
Between classes............ 2 | 2888 | 1444 | 38.0 
|__| | : 
USC) a ai eee nae Leer 23 12416 pee ne (2=0.577) 
85 classes. Inhomogeneity not suggested. 
a MCC.—-x SPEC. 
| 
eins Degrees of Weighted Sum Mean | : 
Variance : : S.D 
Freedom of Squares Square 
re . ah Poe % mice a Hn | ; 
Within classes............. 21 | 8665 413 | 20.3 
Between classes............ 2 | 7 ae, ty 1.87 
IRD Se racic ipcenie eks 23 BOF Ye ne ....| (= —2.385) 
| 
8 classes. Inhomogeneity not suggested. 
mw YERKES—z SPEC. 
cays Degrees of Weighted Sum Mean 7 
Variance : : : S:D:. 
Freedom of Squares Square 
Within classes... ...2...5... 21 16946 807 28.4 
Between classes............ 2 2417 1208 34.8 
RISE RAS aosntciy cs adler are 23 RGROGE. Alcoa sc surie sacait (s=0. 199) 

















8» classes. Inhomogeneity not suggested. 
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TABLE I[]—Continued 


aw SPROUL—7 SPEC. 



































| 
Degrees of | Weighted Sum Mean | . 
/ariance : | ; s S.D. 
Variance Freedom of Squares Square | 
Within: classes. 3. .0..5..... 21 | 12584 599 24.5 
Between classes........... 2 | 4758 2379 48.8 
a — —_——|_____ 
rr | ) 
MOUS (od sm.cstaeh boone 23 | Ck! en Cee ae | (s=0.688) 
8» classes. Inhomogeneity suggested. 
aw GREENW.—7 SPEC. 
Soo Degrees of Weighted Sum Mez . 
Variance Saeed ee meses S.D. 
Freedom of Squares Square 
Within classes... 0.2... 50: 21 7886 376 19.4 
Between classes............ 2 1392 696 20.4 
ROCA ex cicticmca em eines 23 OF78~ foe ae ee (z=0. 308) 
8" classes. Inhomogeneity not suggested. 
aw YALE PHOT.—7z SPEC. 
“4 Degrees of Weighted Sum Mean 
Variance : ‘ : S.D. 
Freedom of Squares Square 
Within classes............ 21 g106 434 20.8 
Between classes........... 2 6821 3410 58.4 
OLB ho exe eet een 23 P E808% Pde ceeurcaee (s=1.033) 
8 classes. Inhomogeneity strongly suggested. 
a CAPE PHOT.—7z SPEC. 
—————————————— —wI 
rae re Degrees of W eighted Sum Mean SD. 
Freedom of Squares Square 
Within classes............ 21 8433 402 20.0 
Between classes........... 2 1822 gl! 30.2 
ROU A oso. ccshaeon rain 23 10255 | 2s RST eae Z=0.412) 








8» classes. Inhomogeneity not suggested. 
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The results of the tests for homogeneity, applied to the trigono- 
metric minus spectroscopic parallaxes for individual observatories, 
may therefore be summarized as follows: The hourly mean differ- 
ences, 7 Mt.W. 60—7 Spec., show no significant dependence upon 
right ascension. The hourly mean differences, a Alleg.—z Spec., 
show no significant dependence upon right ascension. The hourly 
mean differences, t McC.—z Spec., show no significant dependence 
TABLE IV 


METHOD OF CLASSIFICATION ‘ 
PROBABILITY OF 
Data r : | EXCEEDING 2 

FROM 25 CLASSES 


2h Classes 4» Classes | 8» Classes 
aw Alleg.—a McC...... 0.358 0.404 0.565 0.13 
ax Mt.W. 60—7 Spec. 176 | 300 0.122 209 
a Alleg.—a Spec.......| 047 113 OST | 44 
mx McC.—7 Spec.... 010 | 165 —2.385 | 40 
w VYerkés—a Spec... ....| = 6130 =— .058 0.199 67 
ax Sproul—vz Spec.......| 372 560 0.688 12 
aw Greenw.—7 Spec.... .| = i055 | .073. | 0.308 1577 
ax Yale Phot.—7z Spec.. . | 357 663° T0332 | 13 
a Cape Ph.—7z Spec.... 0.747 | —0.014 | O.412 0.01 
Pos. 1 per cent point. . 0.719: | O2722. 3] 0.877 0.01 
Pos. 5 per cent point... 500—Csi/|W 0.510 =—*|| 0.622 05 
Neg. 5 per cent point. . — 513 | —o.761 | —1.484 | 95 
Neg. 1 percent point. . . | —0.740 —1.131 | — 2.300 0.99 
upon right ascension; the value s= — 2.385 for eight-hour classes, 


ascribable probably to a fluctuation of sampling, is of the wrong sign 
for significance. The hourly mean differences, + Yerkes—z Spec., 
show no significant dependence upon right ascension. The hourly 
mean differences, t Sproul — 7 Spec., exceed the 5 per cent points in 
the larger classes, and systematic errors depending upon right as- 
cension are therefore suggested, although with no very great certain- 
ty. The errors, if present, are probably of long period (longer than 
six hours). The hourly mean differences,  Greenw.— 7 Spec., show 
no significant dependence upon right ascension. The hourly mean 
differences, 7 Yale Phot.—7 Spec., exceed the 5 per cent point for 
the four-hour classification, and the 1 per cent point for the eight- 
hour classification; systematic errors of long period in the right as- 
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cension are very probably present. The hourly differences, 7 Cape 
Phot.—7 Spec., exceed the 1 per cent point in the two-hour classifi- 
cation, but yield an insignificant z in the larger classes; the value for 
four-hour classes is negative. Systematic errors, of short period in 
the right ascension, are strongly suggested. If the Sproul, Yale pho- 
tographic, and Cape photographic results are removed from Table 
III, the averages of the remaining entries in the z columns are 0.066, 
0.176, and —o.102, respectively. The standard deviation of a single 
entry in one of the z columns is 0.3 (roughly correct; in the last col- 
umn it is an underestimate), the standard errors of the means are of 
the order of 0.12, and thus the remaining observatories do not show 
any significant collective’ evidence of systematic errors depending 
upon right ascension. 

The last column of Table IV contains the probability of obtaining 
values of s algebraically larger than the observed values, for the 
two-hour classes. They were computed by taking advantage of the 
fact that for the two-hour classes, the distribution of z is very nearly 
normal and of standard deviation 0.31. In the absence of errors de- 
pending upon right ascension, the probabilities should be distributed 
uniformly and at random between o and 1, and they are so dis- 
tributed if the values for Sproul, Yale photographic, and Cape photo- 
graphic are ignored. The general mean 1s 0.43 +0.12 (m.e.), since 
the expected standard deviation of each value is 1/(12)?; the ob- 
served mean is thus in good agreement with the expected mean, 0.50. 
It will be seen that, except for the case of  Alleg.—a McC., the re- 
sults of the present tests are in substantial agreement with those of 
the harmonic analysis. 


THE POSSIBILITY OF RIGHT-ASCENSION ERRORS IN 
THE SPECTROSCOPIC PARALLAXES 

It is desirable to make certain, if possible, that the cause of the 
systematic errors in the few comparisons between trigonometric and 
spectroscopic parallaxes for which such errors are suggested by the 
preceding investigations is not in the adopted system of spectro- 

8 Sometimes statistical results, individually insignificant, may be significant collec- 
tively. 
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scopic parallaxes. Such systematic errors, depending upon right as- 
cension, would if real be present in all the spectroscopic parallaxes 
together that were used in the comparisons, and can therefore be 
looked for in the material of van Maanen’s Table V, as follows: We 
prepare a table of the mean differences between the trigonometric 
and the spectroscopic parallaxes. In this table each row is assigned 
to one of the eight observatories included in van Maanen’s Table V, 
while the columns correspond to the right ascensions. If the hourly 


TABLE V 


| 
| Weighted Sum | 





vise Degrees of Mean | ; 
Variance : | | : | S.D 
Freedom of Squares | Square 
Within classes............ 42 | 43900 | 1045 32.3 
Between classes............| 5 0227. | 245 55.3 
| | 
Total. 47 | FOUST) Pree; | (3=0.089) 
| 
TABLE Va 
| 
we Degrees of | Weighted Sum | Mean : 
Variance : : | : SD. 
Freedom | of Squares | Square 
| 
cue ey ee aie || ies pede of wakes é - 
Within Classes. ss o0desycicecss 21 44751 | 2131 40.2 
Between classes............ | 2 5376 2688 51.8 
Lo" on ee ee 23 soney. . facies wees (@=0.114) 
| 


differences are combined into four-hour classes, there will be six 
columns; if into eight-hour classes, three columns. Means of the col- 
umns are taken, and inhomogeneity is looked for as usual between 
the columns. Weights are used proportional (as usual) to the num- 
ber of stars for each entry. It did not appear worth while to take 
into account the differing accuracies of the various observatories, be- 
cause the standard deviation of the difference between a trigono- 
metric and a spectroscopic parallax does not really differ much for 
different observatories (the ratio of the largest to the smallest esti- ’ 
mate of such a standard deviation being but 1.5 for the eight ob- 
servatories considered); because, moreover, the results of tests of 
homogeneity are never sensitive to small changes in the weights, and 
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because there is no rigorous procedure by which such different vari- 
ances can be allowed for (this difficulty would not, by itself, be very 
serious). 

No significant results are obtained. If we take four-hour classes, 
the results are as in Table V, while the positive 5 per cent point is 
0.45. With eight-hour classes the results are likewise insignificant 
(see Table V a) since the 5 per cent point for this manner of classifi- 
cation is 0.62. We conclude that there is no evidence for the existence 
of systematic errors, depending upon right ascension, in the common 
system of spectroscopic parallaxes used in the comparisons. 

HARVARD COLLEGE OBSERVATORY 
August 30, 1934 








THE SPECTRUM OF P CYGNI1 
By OTTO STRUVE 


ABSTRACT 


Measurements are given of 158 absorption lines, of which 24 are unidentified, and of 
63 emission lines, of which 15 are unidentified. 

The H lines are found to be relatively sharp and strong, their contours suggesting 
the absence of Stark broadening and the presence of turbulence broadening. The pres- 
sure in the H atmosphere must be very low. 

The intensities of the absorption lines suggest a spectral class similar to B1, but the 
lines of C 1 and Mg 11 are much too faint and those of S m1 are too strong. The absorp- 
tion spectrum has certain points of resemblance to 17 Leporis and to the B5 component 
of 8 Lyrae. 

The intensity gradient of the Het lines is steep, suggesting active turbulence. The 
high-level transitions in NV 1 and O 11 are weakened relative to the low-level transitions. 
This may be a result of departures from thermal excitation, although a similar but small- 
er effect is observed in ordinary B stars as a result of the Boltzmann law. 

Recombination is probably responsible for the lines of H and He 1, and possibly for 
Nu. For atoms of high ionization potential such as Si 1v, N 11, etc., recombination is 
negligible. 

The weakness of Mg 11 in P Cygniand in other Be stars is attributed to the depletion 
of the star’s continuous spectrum at the Lyman limit of H. Similarly, the weakness of 
C 11 is caused by depletion at the ultra-violet limit of He t. 

The prevalence of low-level transitions in O 11, and especially in N 11, is regarded as 
an argument against recombination in these ions, and in favor of fluorescence. 

The absorption lines are characteristic of a higher level of excitation than the emis- 
sion lines. Thus, $7 Iv (45.0 volts), C m1 (47.7 volts), and N mt (47.4 volts) occur only 
in absorption. The limit is probably at He 1 (54.2 volts). The absorption spectrum is 
roughly Bo—B1, while the emission spectrum is more nearly B2-B3. 

The elements of high ionization potential are found to originate at greater depths in 
the shell than the elements of low ionization potential. This suggests that the density 
of the shell decreases outward slower than 1/K?, or that the dilution constant increases 
more rapidly than R?. 

A strong line at \ 4396 agrees in position with an O11 line. The possibility is dis- 
cussed that this line is excited by Bowen’s mechanism of abnormal fluorescence. In this 
case a line of He rat \ 515.60 agrees closely with O 11 515.50. No other cases of abnormal 
fluorescence have been found. 

The radial velocities indicate that V increases with the intensity of the line and de- 
creases as we pass from elements of low ionization potential to those of high ionization 
potential. This suggests an accelerated outward motion of the shell. 


I. INTRODUCTION 


There have been a number of investigations of the spectrum of 
P Cygni. After the discovery of bright lines by E. C. Pickering,’ 
J. E. Keeler,? Miss Maury,’ and Vogel and Wilsing* published de- 


* Nature, 34, 439, 1886. 2 Astronomy and Astrophysics, 12, 361, 1893. 


3 Harvard Ann., 28, 101, 1897; 76, 31, 1916. 
4 Pub. Potsdam A p. Obs., 12, 13, 1899. 
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scriptions and identifications of its stronger lines. A Belopolsky® 
identified nitrogen and obtained the displacements of a few lines. 
E. B. Frost® and P. W. Merrill’ made accurate measurements of the 
radial velocities of the stronger lines. R. H. Curtiss* also obtained 
accurate radial velocities of numerous lines. C. S. Yii,? B. P. Gerasi- 
movic’’ and J. Dufay" investigated its continuous spectrum and de- 
rived a very low color temperature for P Cygni. C. S. Beals’ recent- 
ly applied the method of Zanstra to the emission lines of H and of 
He and derived a temperature of about 20,000° K for the exciting 
star. C. T. Elvey's published a set of accurate measurements of the 
total absorptions and emissions of the more prominent lines, together 
with their radial velocities. Still more recently W. J. Williams's gave 
a brief summary of his measurements of radial velocities. 

In spite of this abundance of information collected in the course 
of nearly half a century, the problem of P Cygni is far from its solu- 
tion, and a new contribution toward its elucidation is particularly 
desirable at this time. Within the last few years considerable ad- 
vances have been made in the interpretation of the spectra of novae 
and Wolf-Rayet stars on one side, and of ordinary Be stars on the 
other. A new point of view, not available to earlier investigators, 
can therefore be taken. The present work is an attempt to test the 
theory of an expanding nebular shell in the case of P Cygni. 

The fact that P Cygni is itself an old nova (discovered in 1600) 
renders it especially interesting. Furthermore, its absorption lines 
are not only numerous, but fairly sharp, and should therefore yield 
reliable radial velocities. Similarly, its emission lines, though not 
nearly as narrow as its absorption lines, are narrower than the emis- 
sion lines of most recent novae, and are therefore suitable for meas- 
urement and identification. Finally, the character of its spectrum 

SAp. J., 10, 319, 1899. 7 Lick Obs. Bull., 6, 156, 1911; 8, 24, 1913. 

6 Tbid., 35, 286, 1912. 8 Pop. Astr., 22, 133, 1914. 

9 Pub. A.S.P., 39, 118, 1927. 

10 Harvard Bull., Nos. 852, 857, 867, 1927-28. 

"C.R., 194, 1454, 1932; J. des obs., 15, 45, 1932. 

2 Pub. Dom. Ap. Obs. Victoria, 6, 95, 1934. 3 Ap. J., 68, 416, 1928. 


"4 Pub. A.A.S. (44th meeting), p. 377, 1930. For a complete bibliography containing 
forty-three references see Merrill and Burwell, Ap. J., 78, 87, 1933, star No. 338. 
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suggests a stage of excitation which is not usually found in the spec- 
tra of nebulae, novae, or ordinary Be stars. 

The immediate incentive in beginning this study was furnished by 
the similarity of P Cygni, in many respects, to two other stars re- 
cently investigated at the Yerkes Observatory. It had been noticed 
many years ago by Frost® that certain spectroscopic features of P 
Cygni resembled those of 8 Lyrae, and I have recently found con- 
vincing confirmation of this.'’ No less interesting is its similarity to 
17 Leporis,"° a star containing only one emission line of the P Cygni 
type (HB), but possessing a set of absorption lines the behavior of 
which in many respects duplicates that of P Cygni. 


II. OBSERVATIONAL MATERIAL 


Two excellent spectrograms taken on September 12.25, 1931 (by 
Struve and Morgan), and on October 30.10, 1934 (by Henyey), with 
the single-prism Bruce spectrograph attached to the 40-inch refrac- 
tor were measured for the purpose of identifying the spectral lines. 
Both spectrograms are on the high-contrast Eastman Process emul- 
sion, and they record the photographic spectrum from about A 3925 
to \ 4800. Only two strong lines were measured outside this range. 
The results of the measurements are in Table I. The estimated in- 
tensities are on an arbitrary scale, which is not identical with that 
used for 7 Scorpii’? in column 6. 

Table II contains all elements for which a search was made, with 
their ionization potentials and with the numbers of lines actually ob- 
served. The total number of absorption lines measured in P Cygni is 
110, of which to are unidentified. The number of emission lines is 
43, of which 9 are unidentified and of which all but 3 have corre- 
sponding absorption lines on their violet sides. 

Special attention should be called to the unidentified lines, some 
of which are quite strong. They are not listed among the permitted 
or forbidden lines of atoms and ions which may normally be expected 
in P Cygni. 

Plate II shows enlargements of spectrograms of P Cygni (October 
30, 1934) and of 7 Scorpii (July 7, 1931). 

15 Observatory, 57, 205, 1934. 

Ad. J., 76, 85, 1932. 17 Struve and Dunhan,, ibid., 77, 321, 1933. 
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TABLE I 


SPECTRAL LINES IN P CYGNI 








| ABSORPTION 












































| 
| EMISSION INTENSITY | 
leery rama — Ee ie ss 
r | | | | EMARKS 
Int; | Wek | Into | Vel- TSco | Lab. 
7 | | | 
| HEP? 33-8 
| | | 
3970.05 eel 20 —140 | 30 = 4 8? 
ATOR.FS: icc: 25 184 | 40 | 30 30 
AZAG:, GO 0.00 20 195 | 40 | 6 40 | 
4861.34......| 20 —200 | 30 ae | oie e 
por I 
| Het EP. 24:5 
| 
4026: $3:.. 2 | ==822 I SEG. foxes st I 
3964.73. 10 | 167 20 «| — 33 «I 8 4 
4009 . 27 8 120 2 |— 6 | 4 I 
4026.19 15 153 1 f= & | 30 5 
4120.85 6 114 5 — 3 | 8 3 
4143.77 9 115 2 Pee 10 2 
4168.97 I 83 ° 6 I 
4387.93 10 123 6 b= so | 20 3 
4437-55-- 2 104 [ee } 2 I 
BATT AB 5 5s 15 ro | @ j— & I 20 6 
QTL GO: 6 2 7 aa | 0 j= 7? | 8 3 
AO2T 03 52.5..- 10 —154 0 | 6 fenessc- 4 
_ os ee nl 
| He tt LP. 54.2 
} 
an ae eer ok : 
4685.81 | or | —118 | RN erie 8 |........| Uncertain 
| C mw EP... 243 
| eae ae ee 
| | | | 
A207 2050405. | I | —138 2 | — 38 | 9 | 10+8 | 
Cut LP. 47.7 
Pn ae | es eens 
4187.05 : I eres ° ae | 5 | 10 
4647.40. | I Byes wanes ° sate 9 ite) 
AGSO.36. =." I ° or 5 | fe) BL. Ol 
4051.35. ne oreo ° 5 | 8 Bl. O11 
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TABLE I—Continued 
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| 
ABSORPTION | EMISSION INTENSITY 
A ] a : 5 REMARKS 
| | 
Int. | Vel. | Int. | Vel. T Sco | Lab. 
N11 LP. 29.5 
-- _— 
ZO8 5.085505 4 —113 I + 3 I 6 
3995.00.... 9 110 9 — 453 IO 10 
4035.09 I gI I + 20 3 4 
4041 . 32. I 79 ° I 5 
4043 .54.. I e) 3 3 
4082.28.. I 71 ° I 2 
4227.83.. I 95 O. Whoses ews 2 3 
4236.98... | “O=1 So biees ° 4 6 
4241.80 I 142 ° 4 8 
4426.05?. I 82 fe) ° ° 
4441.99 I 38 OG. Weegezade I 3 
4447 .04 2 122 2 — 14 3 10 
4539.37 ol . ° I 5 
4001.49 4 98 4 + 4 2 8 
4007.17 5 g2 3 + 18 2 7 
4613.88 3 86 2 Broeail 2 6 
4621.40 5 106 4 + 12 2 ” 
4030.55 7 100 9 + 15 5 10 
4043.11 8 — 87 7 + 8 (1) 8 Bl. N II 
Nii IP. 47.4 
| | 
BOO7 2ST 5.506% 2 — 75 r | + 34 | 6 Io Bl. OT 
4379.09 I O. Pace cacves 3 10 
4510.92.. 2 54 fe) 3 6 
4514.80... I 50 o. | 3 7 
4634.16.... I — 46 o. | 2 8 
4040.64... 9 oi a 3 10 Bl. NI 
Ol LP. 34.9 
Ley eee eee 4 — 67 Oo We itd s . 3 10 
4009.90... 2 83 On Beweilws 7 6+4 
AO72.11.. 2 39 fo) 5.5 Dan 6 8 
4075.87. 3 63 6: lowsa cake 8 10 
4078.86 I 51 RE “Whe gssea ahs 3 4 
4085.12 I 72 fe) 3 3 
4119.22. 2 Cr - resect ster 7 8 Bl. He 1 
Pye ey a I O- Treg. sracingus 5 6 
4033/31). | I — 60 Os Neasccpert 5 7 
AISS AS 66% C2 eee O feadannns 4 8 
AIS0:. 70: .. Ro leds anes Oo Perxcerces 6 IO 
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TABLE I—Continued 




































































ABSORPTION EMISSION | INTENSITY 
x cern: Guaieties | - | ans REMARKS 
Int. Vel. Int. | Vel. T Sco | Lab. 
O 1 I.P. 34.9—Continued 
| | | 
4303.92... | o-! [ove fo) Lovet crate 7 5 
4317.60. | = |) es ° 6 8 
AZO: OSs 0-2-0: | 2 | go ° SRR ee a 8 
4345.57 1 | 107 1 | + 20 10 7 
4349.44 5 | 86 fe) 12 8 
4350. 20~. I | 49 ° 8 6 
4306.91. 3 76 o-! + 46 8 a 
4369 . 28? 2 | (4) I (+109) I 4 
4378.40 1 | 74 ° we sawnad ; a 3 
4395.95 5 (144) | 4 ae ae 7 
4414.89 2 71 | I + I! 10 | 10 
4416.97. 2 103 | ° 10 8 
4443.05.. I | saa arch Sess o (| 2 5 
4452.38 I | go | o. | 2 6 
4500.50?... I 142 | o (| I 2 
4590.08.. 2 64 o | 5 9 
4590.19 I 60 oi 5 8 
4638 .86.. 2 63 Oo | 3 6 
4649.15 5 82 | o (| 6 | 10 | BI. Cll 
4061 .65 3 69 | o (| 4 | 9 
4076.25. 2 104 fe) | 3 | 8 
4705 .30.. ies aes ° 3 8 
4710.04?.. 1 | — 74 ° I 5 
| | | | 
Vg tt LP. 15.0 
| 7 | | | 
7S on 2. 2 | —118 | I — 3 | 8 100 | 
| | | | 
| | 
5¢ HE EP: 33.3 
7 | 
4552.61 8 a ae S pers 7 | 9 | 
4567.8 7 103 | @- Bivecases 7 7 | 
4574.75 5 — 92 | o. | 6 | 4 | 
| | | 
52IV LP. 45.0 
4088 . 86. 4 — 60 ° | 10 10 
4116.10 4 =. 40 ° | 10 8 
4212.38 I ' ° 4 3 
4954.14.. Be diem se ote at ° | 4 4 
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TABLE I—Continued 
































ABSORPTION EMISSION INTENSITY 
Xr a et a= ‘ii ee a es REMARKS 
Int Vel. Int Vel. | 7 Sco. | Lab. 
| 
Posal tE (Si Sante Pe 
| | 
4162.64...... ‘ |= 2 ° | 2 10 BLC 
| | | 
ST EP: 346 
| | l 
3928.50 I — 98 ° - 9 | 
3983.76 I 2 fo) I ‘a 
4253.51 4 85 ° 12 10 | BL OW 
4285.00 2 118 fe) : 8 | 
4361.57 I Q2 fe) 3 mo} 
4304.77. I — 108 fe) 2 5 
Call I.P. 11.8 
| | | 
3933.68 | I — 133 ° | | 200 Stellar 
3933.68 A 9 |+ 14 o (| 200 Interstellar 
3908 . 49 oad 3 | — 86 ° 1 150 Doubtful 
stellar 
| Fell LP. 16.5 
| 
| | | 
AQSS 17? >... : I — 65 | ° o | 10 





III. COMPARISON OF P CYGNI WITH NORMAL B STARS 


The spectrum of P Cygni is listed in the Henry Draper Catalogue 
as Bip. Its absorption lines may therefore be compared with those 
of three normal stars, 7 Scorpii (Bo), 6 Canis Majoris (Br), and y 
Pegasi (B2). P Cygni exhibits a number of interesting features: 

1. The H lines are extremely sharp; in fact, they are much sharper 
than in any normal B-type giant. Curtiss and Gerasimovi¢ have re- 
ported diffuse absorption borders on the red sides of the emission 
lines. In the case of Hy this must be due to the group of strong O 11 
lines AX 4345-4350; at He a similar effect is caused by O 11 3973; at 
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H6 a faint indication of absorption may be caused by O It 4105, al- 
though I have not measured this line. The appearance of the 7 lines 
suggests complete absence of Stark broadening. 

2. The He 1 lines are exceptionally strong, surpassing in this re- 
spect ordinary Br stars. 

3. He 11 4686 is almost completely absent. This is in agreement 


with the assignment of P Cygni to class Bt. 


TABLE Ia 


UNIDENTIFIED LINES 








| 
| 
| 
| | MP Cyg \(P Cyg) 
A(Star) Int. Int | Int. 
| 


Abs Em 
| 3957.20 2 | oO 
| 4003.42 3 4004 .62 2 
4020.44 2 4022.12 3 
4039. 24 oe. 4037.13 2 4038.74 |} 1 
4164.07 bs 
4352.47 I 4352.07 | 2 
| 4372.66 I | - 
| 4382 I 4383 .07 ve 
4419.62 3 || 4417-42 7 1419.24 a. 
4431.03 2 4420.04 n 4430.91 Lanea 
4501.14 I 
| 4512.30" o-I 
|} 4010.45 I Vee g es 
* Uncertain. 


4. Si iv is fairly strong, which would place P Cygni between Og 
and Ba. 

5. Or and NV 1 are both strong, but exhibit peculiar intensities 
which are discussed separately in section v. 

6. The great strength of the S 1 lines is unusual. 

7. Equally unusual is the weakness of C 11, which is very strong 
in all three normal stars, and of MWg 1, which is fairly strong in + 
Scorpii and in 8 Canis Majoris, and very strong in y Pegasi. 

It is fairly obvious that the absorption spectrum of P Cygni does 
not correspond to that of a normal B star. On the other hand, there 
is a definite resemblance to the absorption spectrum of the B5 com- 
ponent of 8 Lyrae.'* The latter corresponds to a lower stage of ex- 
citation, and the absorption lines of Ou, Vu, Sm, N 1m, etc., are 


8H. Pillans, 7bid., 80, 51, 1934; Struve, Observatory, 57, 265, 1934 
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TABLE II 


LIST OF ELEMENTS \ 








No. oF Lines | 
ELEMENT LP Pe aa | REMARKS 
Abs. Em. | 
eee 4 4 | Very strong 
i: | 24.5 12 11 | Very strong; all lines present except AA 3935 (1), 
| 4023 (1), and forbidden transitions 
Het 5452 I? © | A 4686 is extremely faint or absent 
Cir. | 24.3 I 1 | \ 4267 is surprisingly weak 
Cit. 47.7 4 o | Strongest lines are very weak 
Nt. | 290.5 19 10 | Strong } 
Ni | 47.4 6 2 | All strong lines are present 
Ox. | 34-9] 34 4 Many lines in absorption 
Om. | 54-9 ° ° d 3961.59 is absent 
New | 40.8 ° | o | Absent 
Mgt | 15.0 I 1 | A 4481 is surprisingly weak 
Al ill | 28.3 ° o | Absent 
Si It. | 16.3 ° o | Absent 
Si mr | 333 3 1 | Strong in absorption; extremely weak or absent 
, | | _ in emission | : 
STIV | 45.0 4 o | Fairly strong in absorption 
Pea. | 19.8 ° o | A weak line suspected at \ 4091 is probably not 
P 11; \ 4475 is absent 
Pu 30.3 fe) © | AA 4080, 4222, 4247 are absent } 
Sin. 23.3 I fe) d 4162 is extremely weak in absorption 
Sul | 34.9 6 o | Surprisingly strong 
Cli... 23.4 ° o | Probably absent; a weak line was suspected at 
| 4572 but this may not be real, or may be due 
| to Cail 
Ki.. ar29 ° o | Absent 
Call.... 11.8 3 o | Stellar K and H (the latter doubtful); inter- 
| | stellar K 
Cait :.. | 51.0 fe) o | Absent 
Oc... 12.8 ° o | Absent 
Ti. 13.6 ° o | Absent j 
verrs|,.. cies ° o | Strongest two lines in this range are 4005.71 (60) 
| and 4023.39 (50). These were suspected to be 
| | identical with the stellar lines 4003.42 (abs.); 
| 4004.62 (em.); 4020.44 (abs.); 4022.12 (em.). 
| However, the normal velocities of —100 
| | km/sec. for absorption lines and o km/sec. for 
emission lines leave a discrepancy of 1 A. 
Weaker laboratory lines of V m1, AA 3952 (4); 
4036 (4); 4183 (35); and 4202 (35) are missing 
in P Cygni 
Crit. is/<P6.6 ° © | 4559 is absent 
Mnu | 15.7 fo) o | Absent 
Feu | 16.5 I © | \ 4233 is probably present 
Nin... 18.2 ° ° Absent ? 
SPAT. os:scekeh DESO ° o | Absent 
PA 5 ..| 14.0 ° o | Absent 
Unidentified .|...... 10 9 | 
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missing. However, the character of the H lines, the exceptional 
strength of the He1 lines, and the weakness of Mg ti and C 1 are 
similar features. 

Some of the strongest unidentified emission lines in P Cygni have 
been observed by Merrill’? in BD+11°4673. These are AA 4383, 
4419, and 4431. 

The gradient effect?’ in P Cygni is very pronounced, thus resem- 
bling such peculiar stars as 17 








P Cyg 

Leporis and the Bs compo- 167 mone 

nent of 8 Lyrae. This is strik- yn te. ear 

° ‘ * a) wee 

ingly illustrated by the inten- | - 

sities of the Het lines in +r i x 

Scorpil, y Pegasi,* and P io! a —_ 

Cygni. ‘The estimates for + / 

“ee 2s ; 8 - 3 

Scorpii and y Pegasi were | p 

made on Mount Wilson coudé = © ae | 

spectrograms, and their scale 4' | 
| 


is approximately twice as ex- 
tended as that used for P | 
| 


Cygni. Theintensities plotted o— 1 11.tttititiii it | 
° 2 4 6 8 17o)0C UC 12s SO16 


N 


in Figure 1 have been reduced 
Fic. 1.—He1 in P Cygni (ordinate), 7 Scorpii 


to the same scale. It is clear, : a 
(abscissa, @) and y Pegasi (abscissa, o). 


both from Figure 1 and from 
Plate II, that the strong He 1 lines predominate in P Cygni, while the 
weaker lines are much more conspicuous in 7 Scorpii and y Pegasi. 

While the gradient has not yet been measured, I estimate it to be 
approximately similar to that of 8 Lyrae (Bs5) and of 17 Leporis. 
Consequently, the total absorption is proportional to the number of 
atoms, for weak and moderately strong lines, while for very strong 
lines the total absorption is almost independent of the number of 
atoms. 

In a recent paper on the gradient effect*® it was pointed out that 
turbulence could explain the observed phenomena. It is probable 
that the same explanation holds for the case of P Cygni. The ab- 
sorption lines, as well as the emission lines, are believed to originate 

9 Ap. J., 69, 330, 1920. 

20 Struve and Elvey, 7bid., 79, 409, 1934. 21 Tbid., 74, 225, 1931. 
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in an expanding nebula. ‘The emission lines are appreciably broad- 
ened. For He 1 4472 I find AX=2.8 A. The contour of the bright 
line is fairly flat at the top, but the edges are not perfectly sharp, as 
would be the case if there were no dispersion in the velocities of ex- 
pansion (or turbulence). In fact, Beals’? has shown that the amount 
of dispersion can be determined from the contours. 

The shape of the stronger absorption lines in P Cygni confirms 
this view. They are appreciably widened and resemble contours pro- 


duced by turbulence broadening.” 


IV. THE INTENSITIES OF THE ABSORPTION LINES 


The absorption lines of O 11 and NV 11 are sufficiently numerous for 
a detailed discussion. Table III gives the more important lines of 
these elements, arranged according to multiplet. The corresponding 
atomic transitions are shown in Figures 2 and 3, which are taken 
partly from Becker and Grotrian’’ and partly from Bacher and 
Goudsmit.?4 The multiplets of both spectra may be roughly subdi- 
vided into two groups: in O11 the group of lower excitation starts 
from configurations 2s’2p’3s (limit *P,,,.) and 2s’2p’3s (limit 'D,) 
and ends in configurations 2s’2p’3p (limit 3P,,.) and 2s’*2p’3p (limit 
'‘D,). The lower states of these transitions connect directly with the 
low metastable states 2s’2p’2p|’P and *D], and, in a few cases, with 


22 In this connection reference may be made to Menzel’s recent criticism of our hy- 
pothesis of turbulence (Pub. A.S.P., 46, 216, 1934). The formulae used were origi- 
nally derived for the case of thermal Doppler effect. Their extension naturally presup- 
poses that the turbulence varies along the radius. This is entirely reasonable, not only 
in expanding nebular shells, like those of P Cygni, 8 Lyrae, and 17 Leporis, but also in 
giant atmospheres like those of « Aurigae and ¢ Aurigae, which have a thickness of the 
order of one astronomical unit. The effect of convection currents (or of turbulent veloc- 
ities which are essentially constant along the radius) has been investigated in connection 
with our work on stellar rotation and found to be inappreciable, at least in those spectro- 
scopic binaries in which a direct test of axial rotation is possible. The broadened lines 
in stars having large intensity gradients are not, however, produced by rotation. The 
contours of lines in 17 Leporis, e Aurigae, etc., cannot be reconciled with rotational or 
convectional broadening. They are in good agreement with the theory of turbulence, 
in the limited sense specified above. 

23 Ergebnisse der exakten Wissenschaften, 7, 8, 1928. 

24 Atomic Energy States, New York, 1932. For O11 see Russell, Phys. Rev., 31, 27, 
1928. The multiplet designations in Table III are from C. FE. Moore, A Multiplet Table 


of Astrophysical Interest, Princeton, 1933. 
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the ground level of the O Irion, 2s?2p’2p'S. The upper levels of these 
transitions do not directly connect with the low 2p-levels. In the 
second group the transitions, corresponding generally to higher-ex- 
citation energies, start from those levels which were upper levels in 
the preceding group and end in transitions 2s*2p’3d (limits *P,,. and 
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'T),). Finally there are a few lines, also included in the second group, 
which start from configurations 2s*2p’3d and end in configurations 
2s°2p’4f. The average excitation potentials of the two groups are 23 
volts for the lower and 26 volts for the upper state, in the first group; 
and 26 volts for the lower and 29 volts for the upper state, in the 
second group. 

A similar division has also been made in the case of V 11. Here the 
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excitation potentials of the first group are approximately 18 and 21 
volts, while those of the second group are 21 and 23 volts. 

In Figure 4 I have plotted the estimated absorption intensities in 
the spectrum of P Cygni against those of r Scorpii, 8 Canis Majoris, 
and y Pegasi. The first group of lines is represented by filled circles, 
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Fic. 3.—Grotrian diagram for N11 


while the second group is represented by open circles. A few un- 
classified lines of V 11 have arbitrarily been added to the second 
group. The two O u lines, \ 4469 and X 4396, have been omitted for 
reasons which are explained in section vii. Several other lines which 
are complicated by blends have been omitted. These are shown in 
parentheses in Table III. 

In all four sections of Figure 4 the filled circles are systematically 
above the open circles, indicating that for both Om and Nu the 
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upper transitions are weakened relative to normal B stars. This may 
also be confirmed by an inspection of Plate II. 

The suggestion might be made that this is a normal effect of ex- 
citation. Indeed, I have found some time ago’ that the O 1 lines do 
not all come to a maximum at the same place in the spectral se- 


Nu, r Sco Nu, y Peg 


P Cyg 
NO 


Intensity: 


| 

| 

o 2 & & & to 2 Oo 2 4. 6 [St Ace 12) 14 
O11, r Sco Ou, B CMa 


I'tc. 4.—Intensities of high-level transitions (0) and of low-level transitions (@) 


quence, but that those of higher excitation are, in accordance with 
theory, slightly stronger in the earlier spectral subdivisions. 

In order to test this effect, I have compared * in Table IV several 
of the most conspicuous lines of both groups in P Cygni, 10 Lacertae 
(Og), and Herculis (B3). Of the two comparison stars, the first cor- 
responds to a higher general level of excitation than P Cygni and the 
second to a lower level. The relative strengthening of the high-level 
lines, in passing from class B3 to class Og, is obvious. 


*5 Ap. J., 78, 73, 1933- 
* The intensities of lines in stars other than 7 Sco are taken from ibid., 74, 225, 1931; 
those of r Sco are from ibid., 77, 321, 1933. 
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TABLE III 






































| INTENSITY | 
WaAvE-LENGTH MULTIPLET ; | \ | REMARKS 
P Cyg : . 
| eg Lab. | T Sco B CMa | 
Low Transitions of O11 
| — ’ 
3973: - 2P—er 4 10 3 5 
4415. z7P—e?D° 2 10 10 5 
ty. (2) 8 10 4 
52. I 6 2 2 
4317 y‘P—e'’Pe 2 8 6 4 
20. 3 8 7 5 
45. (1) 7 10 3 
49. 5 8 12 4 
67 3 7 8 4 
4640. y*P—eiD° | 5 10 6 5 Biel 
in + Sco 
38. 2 6 3 4 
61 3 9 4 3 
76 2 8 3 3 | 
4351. y2D—f?7D°_ | I 6 8 3 
47 (0) 5 8 2 
4591. y7D—-e’F? | 2 9 5 4 
96 I 8 5 2 | 
High Transitions of Om 
— | | | 
4396. eir~xD (5s) | 7 | 3 3 | | 
69. (2) 4 | I I | 
| | 
4705. ei —7F I | 8 | 3 3. CO 
| | 
46099. e?D°—z4D fo) 7 2 | f 
4153... e4Pe—x4P I 7 5 4 
I 6 5 3 | 
4119. e#Pe—ziD (2) | 8 7 4 | 
4076 e4D°—z‘F 3 | IO 8 6 | 
72. 2 8 6 5 
69. 2 6 ¥ 7 
> | ? 
78 I 4 3 3 } 
85. I 3 3 3 | 
| 
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Considering the strength of some O 11 lines in P Cygni, it is some- 
what difficult to assign to it a class later than Br. On the other 
hand, the relative intensities of the high-level and the low-level lines 
would place it nearer to B2, or even to B3. Such a late spectral class 
is not possible because of the strength of the S 111 lines and the pres- 
ence of N 11, Sz Iv, etc. 

In the case of NV 11 the effect is even more striking. No possible ex- 
citation can account for the enormous strengthening of low-level 


TABLE IV 
OF 
INTENSITY 
WaAveE-LENGTH - ; 
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| P Cyg(Bip) | to Lac (Oo) t Here (B3) 
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Low Transitions 
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4415 | 2 fo) 2 
4317 2 ie) 3 
BRE 6.5 0:5 42:5 | 3 I fe) 
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| High Transitions 

| | 
4076 ’ 3 4 I 
Co)’ | 2 | 2 I 
BODO .....%- 2 4 2 
4795 I I fe) 
4099. . oO I ° 





lines of VN 1rin P Cygni. Table V illustrates this. There is little dif- 
ference between 6 Canis Majoris (Br) and 67 Ophiuchi (Bs5), and 
the enormous enhancement of the low-level lines in P Cygni must 
be due to causes other than thermal excitation. 


V. COMPARISON OF EMISSION LINES WITH ABSORPTION LINES 


It is generally recognized that the most effective mechanism in 
the production of bright lines in stellar spectra is photo-electric ion- 
ization with subsequent recombination. This mechanism has been 
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analyzed by Zanstra,?? Menzel,’* Bowen,?? Miss Payne,*° Beals,** 
Woolley,*? and others. Beals has applied Zanstra’s theory to the 


TABLE V 

















Nil 
| INTENSITY 
Wave-LencTu 
| P Cyg | B CMa (Br) | 67 Oph (Bsp) 
| | 
| Low Transitions 
3995. | 9 4 | 3 
4031. 7 3 | 4 
4613 .. 3 2 | I 
4021 5 I | 2 
| High Transitions (and Unclassified) 
0 SSR Ch RE Ie o-! | I ° 
4448 2 2 2 
4035 I 2 I 
4242 I 3 2 


emission lines of H and of He 1in P Cygni and has found a tempera- 
ture of 19,300° K in the former case and 26,000° K in the latter case. 


27 Thid., 68; 50, 1927;-25. f. AD., 21,1031. 

28 Op. cit., 38, 295, 1926. 29 Ap. J., 67, 14, 1928; 81, 1, 1935. 

39 M.N., 92, 368, 1932. Miss Payne’s criticism of my rotational hypothesis of bright 
lines is of no consequence to this paper. There is no disagreement in so far as recombina- 
tion and fluorescence is concerned. In fact, I have mentioned these mechanisms in 
several papers and have accumulated observational evidence in favor of recombination 
as the principal source of bright H/ lines in Be stars. There are five major questions to 
be answered in connection with stellar emission lines: (1) Where do these lines originate? 
(2) How do the atoms get there? (3) By what mechanism are they excited to radiation? 
(4) Why are some bright lines broad and flat-topped while others are narrow? (5) Why 
do many bright lines vary in intensity and contour? There is no discordance with re- 
gard to question 1. I have attempted to answer questions 2 and 4 by the rotational 
hypothesis for ordinary Be stars, and by the expanding shell hypothesis for stars of the P 
Cygni type. Question 3 is satisfactorily answered by the work of Rosseland, Miss Payne, 
and others. Question 5 has not yet been adequately answered, although McLaughlin 
and Gerasimovié have made interesting attempts in this direction. 


* Pub. Dom. Ap. Obs. Victoria, 4, 271, 1930. 3? M.N., 94, 631, 1934. 
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Since the emission lines of eI are very strong, it is obvious that 
there must be a large amount of energy in the star’s spectrum toward 
the violet from \ 504. This agrees with the temperature derived by 
Beals. The energy-curve corresponding to 20,000° K and a few of 
the series limits are shown in Figure 5. 

N iris only partly in emission, indicating that at \ 418 the energy- 
curve is somewhat depleted. Sz 11 and O 1 appear almost exclusive- 





° 400 800 I 200 1600 2000 2400 A 


Fic. 5.—Black-body curve for 7’ =20,000° K 


ly in absorption, which might indicate an insufficient amount of 
radiation at \ 370 and 353. Elements of higher ionization poten- 
tial, such as S 11, V 11, Sz iv, and C 1, are seen only in absorption, 
while O 11 is completely absent. 

The emission lines agree fairly well with the recombination hy- 
pothesis, but there are a few anomalies. The absence of Sz 11 and the 
weakness of Mg tare interesting. The first, though its ionization po- 
tential of 16.3 volts is slightly above that of 7, would not normally 
be expected in a B2 spectrum, since in ordinary stars the relatively 
low ionization potentials of Sz 111 and Si tv will lift a large number of 
atoms into higher stages of ionization, thus preventing a recombina- 
tion spectrum of Sz 11. However, we have just pointed out that there 
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is probably no recombination spectrum of S7 111, so that the majority 
of the Si atoms must be in the ground level of the Sz m1 ion. Even 
more surprising is the weakness of Mg 1. In normal B2 stars, \ 4481 
is a fairly strong line, and it usually persists into the O’s. Now, the 
ionization potential of Mg 111 is 80 volts, so that the majority of Mg 
atoms must be in the ground level of Mg 111, and should thus favor 
recombination processes forming the line \ 4481 in emission. 

The most obvious explanation of the weakness of Mg 1, a phe- 
nomenon that seems to be prevalent in Be stars, is that the star’s 
continuous spectrum in the neighborhood of \ 824 is still sufficiently 
depleted of radiation by the strong photo-electric ionization at the 
limit of the Lyman series of 7, at \ g11. This would not be surpris- 
ing, since the mechanism of Zanstra, which is justified by its results, 
demands that a large fraction of the entire energy on the violet side 
of \ g11 is utilized in producing the emission lines of 7. 

A similar but even more striking anomaly exists in the case of 
Cu, which should normally be strong in class B2. Here the ioniza- 
tion potential, 24.3 volts, is almost identical with that of He 1, 24.5 
volts, and the limits, at \ 508 and X 504, are so close to one another 
that the recombination spectrum of C 11 will doubtless be suppressed 
by continuous absorption at the ultra-violet He 1 limit. 

Much less satisfactory is the behavior of V 1. Here the ioniza- 
tion potential of 29.5 volts corresponds to \ 418. The amount of 
radiation available for ionization at this wave-length is about one- 
seventh of that available at the He 1 limit and about one two-hun- 
dredth of that available at the 7-limit. At the limit of On, A 353, 
where no appreciable recombination occurs, the energy is one- 
eightieth of that at \ 504, He 1. We are therefore uncertain whether 
recombination does or does not occur in N It. 

However, the observations show clearly that the emission is lim- 
ited to low-level lines, strong high-level lines being devoid of emis- 
sion. This is distinctly unfavorable to the hypothesis of recombina- 
tion: in the laboratory, recombination has invariably increased 
the intensities of high-level lines, and Mohler and Beckner* have 
even succeeded in observing new high-level lines which are not seen 

33 J. of Research, Bureau of Standards, 2, 489, 1929. See also Ruark and Urey, Atoms, 
Molecules and Quanta, p. 470, 1930. 








86 OTTO STRUVE 


in ordinary laboratory sources. The anomaly of the NV 11 emission 
lines is probably related to that of the absorption lines discussed in 
the preceding section. 

The absorption lines of P Cygni are characteristic of a much higher 
degree of excitation than are the emission lines. The limit is prob- 
ably in the vicinity of 50 volts: Sz Iv, ionization potential 45.0 volts, 
and N 111, ionization potential 47.4 volts, are present; C 111, ioniza- 
tion potential 47.7 volts, is very weak. 

It is quite clear that while the absorption lines of P Cygni corre- 
spond to a class intermediate between Bo and B1, the emission lines 
are more nearly at B2 or B3. It might be thought that the appear- 
ance of emission is merely a function of the intensity of the absorp- 
tion lines, but that is not true. Thus, some faint lines of He 1 and of 
N it are accompanied by emission, while some of the stronger lines 
of O 11 are devoid of it. 

The radial velocities discussed in section viii leave no doubt that 
the absorption lines originate in the expanding shell, and in this re- 
spect there is no great difference between absorption lines which are 
accompanied by emission lines and those which are not. 

There is at present no adequate theory of the origin of absorption 
lines in expanding nebular shells. It is obvious, however, that in a 
general way we are dealing with the same process of absorption and 
re-emission (or scattering) which gives rise to the absorption lines 
in a reversing layer. 

Ignoring for the moment all anomalies arising through departures 
from thermodynamic equilibrium, such as were mentioned in the 
preceding paragraphs, we should expect that for a thin shell of radius 
R around a star of radius 7, the ratio of absorption intensity to emis- 
sion intensity would be roughly proportional to the ratio of the vol- 
ume of gas projected upon the star’s disk to that of the entire shell 
visible to the observer. This gives 


2 


rr 
= = 
A NR. 
a ey: 
y yr? 
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Accordingly, if for two equal absorption lines the emission lines are 
not the same, we conclude that they must originate in layers of differ- 
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ent R. The function of (NV, o), depending upon the number of atoms 
and the absorption coefficient, takes care of the gradient effect, etc. 
The formula, of course, fails if R=r, for then we should use the cus- 
tomary treatment of a scattering reversing layer. The factor f 
(V, a) also shows that it is not permissible to compare unequal ab- 
sorption lines. 

Since we have shown that equal absorption lines may have widely 
different emission lines, we conclude that, in general, the former 
originate in other shells than do the latter. This agrees with the re- 
sult that the level of excitation is much higher for absorption than 
for emission. 

This effect is seriously complicated by departures from thermo- 
dynamic equilibrium. These must be made responsible for the enor- 
mous difference in the emission and absorption ratios of He 1 lines 
belonging to different series. Thus, Schwede** has found, in agree- 
ment with Elvey:"s 
He 4388: total absorption... 0.62A He 4388: total emission.... 0.16A 

4472: total absorption... 0.69 4472: total emission.... 1.56 
It is difficult to explain this result, unless the processes of absorption 
and emission are far from being balanced. 

Eddington* has shown that the state of ionization in a nebula is 
equal to that of a reversing layer having a density of p6, where p is 
the density of the nebula and 6 the dilution factor 4R?/r*. Since, for 
the absorption lines of highest excitation, A/F is smaller than for 
lines of low excitation, we conclude that R(high exc.) < R(low exc.), 
and that, consequently, the factor pé decreases as we approach the 
star. Since 6~k’, we find that p must increase inward more slowly 
than 1/k?. According to Eddington,*® in an isothermal nebula 
p~1/R?’. An alternative explanation would be that 6 increases out- 
wardly more rapidly than R?. 

VI. NORMAL EXCITATION OF SPECTRAL LINES IN P CYGNI 

If the absorption lines of P Cygni originate in an expanding shell, 
it may not be permissible to apply to them the ordinary theory of 
stellar absorption lines. Especially does it seem dangerous to apply 
the Boltzmann law for the derivation of the numbers of atoms in 


34 Ap. J., 77, 348, 1933. 
35 Internal Constitution of the Stars (German ed.), p. 478, 1928. 3° Ibid., page 489. 
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various quantum states. The problem of radiative equilibrium in a 
nebula has been treated by Rosseland’’ and applied by him to the 
case of three levels of increasing energy, 1, 2,3. Ifa; and a, are the 
transition probabilities, the condition of equilibrium is 


LaKiXe=O; ix — Ddiz ‘ 
k k 


The solution of these equations gives the numbers of atoms in the 
various states, 

Xi =Amak™ ’ 
where X,, is a constant and a*” is the subdeterminant of the deter- 
minant d;,. For three states we have 


a = Az,(dj: +432) +4230; , 
a" = A;,(d3, +432) +41;032 , 
a3" = d,3;(da +423) +1242; « 


Rosseland applies these formulae to the case where the second level 
is metastable, i.e., a,,=a,,=0, and finds that in this case the number 
of atoms in the metastable level is similar to that given by Boltz- 
mann’s law. 

To represent more closely the conditions prevailing in the lower 
levels of the two groups of transitions of O 11 and of NV 11, we shall as- 
sume that a,,=d,,=0. In that case 


A= 032, A7=Ay2032 , A= 44202; , 


and, accordingly, 
x, a" an 
2 a tie’ 
x; _ ast a Qy2423 


t, G" Gy0y2° 
Now, according to Einstein, for emission a;;=a(1+plvi.|) and for 
absorption Aig = Aix (viz) ./W; 2 
p(vi.) is the energy density of the radiation divided by 8rhv3/c’. 
Accordingly, 


p(vix) =. plviz) , 


37 Astrophysik auf atomtheoretischer Grundlage, p. 224, 1931. 
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where p(v;,.) is the energy density given by Planck’s law. Substitut- 
ing this, we find 


Xz @2 p(v;2) I 


xy Wy 1+ (v2) 





x; @; p(V42)p(V23) I 


rt 5 plval)(t+5 elva) 


xX, 


For small values of p(v) the number of atoms in the second state is 
proportional to 1/6 and the number in the third state is proportional 
to 1/6*. It is obvious that if 6 is large enough, there must be a mark- 
edly smaller number of atoms in state 3 than would be given by 
Boltzmann’s law. Furthermore, it would appear that for suitable 
values of 6, «; is smaller if a,,=a,,;=0 than if a,,=4a,,; #0. 

In a nebula with large 6 we should therefore expect that aside 
from recombination lines and resonance lines (which, for most atoms, 
are not in the observable range), only those lines would be strong 
which have lower levels connecting directly with the ground level or 
with one of the low metastable levels.** This condition is fulfilled for 
Nu and probably also for O11. It may hold in the case of Si 1n, 
whose lower levels are reached directly from the metastable *P-state, 
but it is violated in the case of Sz 1v, whose lines \ 4089 and A 4116 
originate in the 2s’S-state, which does not directly connect with the 
ground level, 1s’S, but can only be reached from it by two or more 
jumps. It is therefore probable that, at least for Sz Iv, 6 is not very 
large. A numerical evaluation of the quantities involved is not, at 
present, possible. We can only conclude that slight departures from 
thermal excitation are present, but that for Si 1v conditions resemble 
those prevalent in ordinary stars. 


VII. ABNORMAL FLUORESCENCE 
The spectrum of P Cygni contains a peculiarly strong line, present 
both in absorption and in emission, at \ 4396. This line has been 


38 We ignore transition probabilities from the ground level to the upper states of the 
high-level lines. Although this is not rigorously correct, it seems probable that these 
transitions will have little effect in lifting atoms into a state from which emission of the 
high-level lines is possible. In the first place, the continuous spectrum in the required 
wave-lengths (A 430+) is relatively weak; in the second place, their coefficients a are 
probably small. 
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tentatively identified as O11 4395.95, and all available measure- 
ments prove that the agreement in wave-length is excellent. How- 
ever, the violet absorption line is distinctly too strong for its labora- 
tory intensity and the emission line, if really due to O 11, is by far the 
strongest line of this atom. 

Two possibilities suggest themselves: either (1) the agreement in 
wave-length is a mere coincidence, the line being really unidentified, 
or (2) the O 11 atom has an abnormal tendency to enhance this par- 
ticular transition. 

It is difficult either to support or to disprove the first possibility. 
There are four other unidentified lines in the region AA 4300-4500, 
and \ 4396 may belong to this group. On the other hand, a chance 
coincidence of an unidentified line with O 11 4395.95 is improbable. 

The second possibility merits a careful scrutiny. Bowen’? has re- 
cently found a powerful secondary mechanism of excitation which 
depends upon the close agreement in wave-length of a strong ultra- 
violet emission line (itself produced by the primary mechanism of 
recombination) with an absorption line of another atom. The ques- 
tion in our case is: Is there such an anomalous mechanism at work 
in Ou which would lift an abnormal number of atoms into either the 
upper or the lower level of the line \ 4396? 

The line Om 4396 is a member of the multiplet 2s’2p?3d?D — 
2s°2p*3p’D. The upper terms are v= 48,618.4 and v= 48,566.4. The 
lower terms are v=71,498.9 and v=71,308.2. The first upper level 
leads to the lines \ 4369.28 (4) and \ 4406.02 (1), while the second 
gives \ 4395.95 (7) and A 4359.38 (1). The two upper levels do not 
connect with the ground level, which is a 4S, but they connect with 
the very low metastable level 2s’2p’2p’P, viz., v= 242,555.5 and v= 
242,560.0. These lead to the two following close pairs: \ 515.62 and 
515.49. Bowen»? has actually measured \ 515.62 (2) and X 515.47 
(2). The first of these feeds \ 4367 and the second feeds \ 4396. 

Bowen has pointed out that only H and He could reasonably be 
expected to produce abnormal fluorescence. H and He 11 have no 
lines near \ 515. But there is a conspicuous He 1 line, 1'S—5'P (see 
Grotrian diagram, Fig. 6), which gives, from the term values, 
515.07 and which has been measured by Lyman’ at \ 515.65 and 


39 Phys. Rev., 29, 243, 1927. Ap. J., 60, 11, 1924. 
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by Hopfield* at \ 515.596. The agreement of this value with the 
two O 11 lines is rather striking, being 0.02 A in the first case and 
o.11 A in the second.*? The second, larger, residual applies to the 
case of O I 4395. 

The visible He 1 lines, especially of the series 2*5 —m'P, are very 
strong in emission. The width of \ 4472 is approximately 2.8 A, and 
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Fic. 6.—Grotrian diagram for He 1 


the contour is fairly flat at the top. Since the broadening is due to 
Doppler effect of expansion,*’ we have 


AAmd , 


4" Tbhid., 72, 136, 1930. 

” Since writing the foregoing, Professor Bowen has called my attention to the meas- 
urements of the O 11 lines by B. Edlén (Nova Acta Regiae Societatis Scientiarum U psali- 
ensis [4th ser.], 9, No. 6, 1933): \ 515.640 and \ 515.498. These values are estimated 
to be accurate to c.or A or c.o2 A. The residuals are therefore reduced to 0.098 A in 
the case of \ 4396 and 0.044 A in the case of \ 4369. There is still a decided preference 
for \ 4369. Professor Bowen also mentions that K. T. Compton and J. C. Boyce (J. of 
Franklin Inst., 205, 497, 1928) found for the He 1 line a value which is 0.004 A higher 
than Hopfield’s more recent result. 


‘3 Beals has found that in the visible region the line widths are proportional to their 
wave-lengths (op. cit., p. 297). 
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and the width of Het 515.60 should be 0.32 A. Allowing a similar 
width for O 11 515.5, we find that both O 1 lines should profit from 
the excessive radiation available for absorption from the bright He 1 
line. It is obvious, however, that unless our computation is in error, 
d 4369 should profit a great deal more than A 4395. Normally \ 4360, 
of laboratory intensity 4, would not be expected to show in P Cygni. 
There is a faint, but definite, absorption line at approximately the 
correct wave-length, and there is also an emission line, but its posi- 
tion is rather markedly displaced toward the red. (A slight displace- 
ment of the absorption line in the same direction is not serious, be- 
cause it blends with a weak bright line at O 11 4367.) I am therefore 
in doubt whether A 4369 is really the O 11 line. Should it not be O 1, 
I should be inclined to abandon also the hypothesis of abnormal 
fluorescence in \ 4396 (laboratory intensity 7) and should prefer to 
regard this line as unidentified. 

Unfortunately, no other multiplets which are easily observable 
arise from the O11 level y= 48,566.4. The multiplet AX 4941, 4943, 
4956 cannot be photographed with the Eastman Process emulsion. 

I have made a search for other coincidences in ultra-violet wave- 
lengths (see Fig. 7). The nearest are: 


Pera sabe ek ie Ot OE i ccd sss RR. BE ce os Sa 
1 oer re i Se eS) ree 537-8 
NE, Sotto’ C4 ae ee Oe 538.3 


The first coincidence may be sufficiently close to produce observable 
effects, but unfortunately the ordinary lines of S 1 are all faint. Not 
one of the unidentified emission lines agrees with classified lines of 
Sr listed by S. B. Ingram.44 The second coincidence is not close 
enough; furthermore, the He line corresponds to a forbidden transi- 
tion and is probably weak in P Cygni. The third case is too remote. 


VIII. THE RADIAL VELOCITIES 
It has been known for a long time that radial velocities obtained 
from different lines of P Cygni are not consistent among themselves. 
Thus, the absolute values of the negative velocities from the H lines 
become larger as we pass from the violet to the red, and for this 
reason the hypothesis has been advanced that V=f(A). For He1 


4 Phys. Rev., 32, 172, 1928. 
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this functional relationship was found to hold when only a few of the 
stronger lines were measured. 

I have made no attempt to obtain accurate velocities from the 
stronger lines, but have measured velocities for as many lines as 
possible. The results are given in Table I and are summarized in 
Table VI, in which there are given for each element the mean veloc- 
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Fic. 7.—Ultra-violet spectra 


ities corresponding to lines of approximately equal absorption in- 
tensity. ‘The elements are arranged in order of their ionization po- 
tentials. 

In spite of the large scatter, there is a fair degree of correlation be- 
tween intensity and velocity; for any given element the strong lines 
show a larger velocity of approach than the faint ones. There can be 
no doubt that this is the true correlation, and that the physically 
absurd correlation with wave-length was merely a result of the gen- 
eral tendency in H and He 1 for the stronger lines to have the longer 
wave-lengths. Exactly the same phenomenon has been observed in 
the B5 component of 6 Lyrae. 
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The most interesting result, however, is the obvious correlation 
between velocity and ionization potential: for two lines of a given 
intensity, the one of lower ionization potential has the larger veloc- 
ity of approach. This effect is not caused by the presence of emission 
lines; thus for C 11 and Mg u1, which are very little affected by emis- 
sion, |V| is much larger than for N 1, which has strong emission. 
Furthermore, the transition from emission to absorption is gradual, 
not sudden. 

TABLE VI 


RADIAL VELOCITIES IN P CYGNI 





| meeenee Average | 
Element | Intensity a Pe No. of Lines 
Velocity 
| (Absorption) | 
| km/sec. 
| te Cee ae a 22 —181 4 
Mgt (15.0) 2 118 I 
hi OY Re) I 138 I 
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He oon AoA 5 dusts ots 2 113 | 2 
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Het crisp tric ech en ia a 8 15 156 | 2 
Nu (29 5). ep Seer I 85 7 
Nu 3 5 | 4 
Niu 5 99 2 
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S7 Tr (33.3) 2 5 g2 I 
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CORE 5) caateiate: Paden aaa es 3 75 | 4 
| ee ae ee ave ans 78 | 3 
OLIV (A550)... 4 5a | 2 
N rt (47.4). bts peat I 48 | 2 
6 a a re ae poaR ates 2 — 64 | 2 
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We have consistently been led to the result that lines of high exci- 
tation originate deeper inside the shell. Consider a single element, e.g., 
Het. The strong lines originate in smaller optical depths within the 
nebula than the fainter lines. The velocities of the strong lines are 
greater than those of the faint lines, suggesting an accelerated mo- 
tion of the nebula outward. Now consider two lines of the same in- 
tensity but of different elements. The one of lower potential, as 
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found in section vi, originates at a greater distance from the star than 
the one of higher potential. Consequently, the former should show 
a larger velocity of approach than the latter. This is in agreement 
with observation. 

There is good evidence in favor of this acceleration, from the 
emission lines. These are by no means all of equal width, but a 
closer study of their behavior must be postponed until actual con- 
tours are available. In passing, however, attention may be called 
to the sharpness of the unidentified emission line \ 4164. 


IX. THE ULTRA-VIOLET AND GREEN REGIONS 

For completeness I have measured several spectrograms on coarse- 
grained, isochromatic plates taken by Mr. J. A. Hynek with the 
prism-spectrograph attached to the 69-inch Perkins reflector. The 
measures are uncertain, especially in the green region, where the dis- 
persion is small and where the comparison lines are not favorably lo- 
cated. The list contains the designations or the laboratory wave- 
lengths of all identified lines and the measured wave-lengths (cor- 
rected to the sun) of the unidentified lines. Attention is called to the 
peculiar intensity of Si 11 3806.56, which is a high-level line. The 
forbidden O 11 lines AX 3726, 3729 are not present, thus confirming 
the conclusion that 6 is relatively small. The O11 multiplet 4941, 
43, 56 (referred to in section vii) is not definitely seen on the coarse- 
grained plates, although the first member may be identical with 
d 4938.01 (unidentified). The list contains 48 absorption lines (of 
which 14 are unidentified) and 20 emission lines (of which 6 are un- 
identified). 

ABSORPTION LINES 

H.; (1) —136; H.; (1) —122; Hx (2) —143; Ha (2) —143; Hoo (2) —138; 
H,, (2) —105; His (3) —134; Hi; (4) —152; Hw (7) —164; Hi; (4) —122; Ai, (4) 
—146; H,; (4) —105; Miz (5) —133; Hus (4) —142; Ho (7) —138; Hy (8) — 1433 
H;=Hs (20 bl. Het) Het: 3587.30 (2) —198; 3613.64 (3) —183; 3634.30 (6) 
— 113; 3705.10 (bl. H); 3819.63 (8) —132; 3867.50 (3) —111; 3871.80 (3) —91; 
5015.68 (7) —180; 5047.74 (1) —176; O I 3727.33 (2) —95; S I: 3632.03 (1) 
— 216?; 3709.37 (2) —725 3717.77 (1) —130; 3778.91 (1) —96; 3860.64 (1) — 47; 
Si 11? 3791.41 (1) —152; St IV? 3773.13 (1) —68. 

Unidentified: 3599.56 (3); 3602.51; 3648.34 (1); 3723.89 (1); 3782.19 (2); 
4938.01 (O 11 4941.12?); 5052 (1); 5071 (1); 5082 (2); 5007 (1); 5125 (2); 5131 
(1); 5130 (2); 5154 (3). 
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EMISSION LINES 

His (1) —17; Hi; (2) +1; His (3) +21; Ais (2) +26; Hi, (2) —14; Hi; (3) 
— 26; Hi2 (3) —17; Au (3) —14; Hiro (4) —14; Hy (6) +0; He 1 3819.63 (6) —10; 
3888.90 (bl. H); 5015.68 (10) —6; Si 111 3806.56 (4) —32. 

Unidentified: 3600.88 (3); 3603.49 (2); 3757-57 (1); 5055 (3); 5127 (2); 5156 
(3). 

X. CONCLUSIONS 

As a result of the preceding discussion the following tentative pic- 
ture of P Cygni suggests itself: 

The nucleus is a star of effective temperature 20,000° K, situated 
at a distance of approximately 1ooo parsecs from the sun. Owing to 
its great distance, the interstellar line K is exceptionally strong. Se- 
lective interstellar absorption reduces the color temperature to ap- 
proximately 6000° K. There is no line spectrum belonging to a 
stationary reversing layer similar to that found in 17 Leporis. ‘The 
absorption lines and the emission lines originate in an expanding 
nebular shell. The velocity of this shell is accelerated outward, be- 
ing approximately 200 km/sec. in the region of formation of the // 
lines and approximately 50 km/sec. in the region where S7 Iv and 
N 101 are formed. 

For a given element the stronger absorption lines originate at 
lesser depths than the weaker lines. For two equal lines of different 
elements the one of higher ionization potential originates at the 
greater depth. 

The mechanism of line excitation deviates from thermal excita- 
tion, and peculiar cases of fluorescence may be present. It is prob- 
able, however, that the dilution constant in P Cygni is much smaller 
than that found for novae or planetary nebulae. 

YERKES OBSERVATORY 
November 1934 

45 Since this paper was written, I have received the November, 1934, issue of the Ob- 
servatory, containing an article by C. S. Beals on P Cygni (p. 319). It is gratifying to find 
that his conclusions in regard to stratification in the shell agree with mine. In this con- 
nection attention may be called to the star 17 Leporis in which the ratio emission /ab- 
sorption is very small and for which I concluded that the ratio: radius of shell/radius of 
star, must be relatively close to 1; this is, of course, the same result as that given by 
O and by other atoms of high ionization potential, in P Cygni (Ap. J., 76, 103, 1932). 
Beals now favors a temperature of 30,000° K. 





